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ABSTRACT 

This paper investigates on the possible systematic difference of Supernovae la (SN la) prop- 
erties related to the age and masses of the progenitors that could introduce a systematic bias 
between low and high redshift SN la's. The relation between the main features of the distri- 
bution of the delay times (DTD) and the masses of the progenitors is illustrated for the single 
(SD) and double degenerate (DD) models. Mixed models, which assume contributions from 
both the SD and DD channels, are also presented and tested versus the observed correlations 
between the SN la rates and the parent galaxy properties. It is shown that these correlations 
can be accounted for with both single channel and mixed models, and that the rate in SO and 
E galaxies may effectively provide clues on the contribution of SD progenitors to late epoch 
explosions. A wide range of masses for the CO WD at the start of accretion is expected in 
almost all galaxy types; only in galaxies of the earliest types the properties of the progenitors 
are expected to be more uniform. For mixed models, late type galaxies should host SD and 
DD explosions in comparable fractions, while in early type galaxies DD explosions should 
largely prevail. Events hosted by star forming galaxies span a wide range of delay times; 
prompt events could dominate only in the presence of a strong star-burst. It is concluded that 
nearby SN la samples should well include the young, massive and hot progenitors that neces- 
sarily dominate at high redshift. 

Key words: stars: evolution - supernovae: general - galaxies: evolution - cosmology: dis- 
tance scale 



1 INTRODUCTION 

Type la Supernovae (SN la) pro vide an important contribu tion of 
iron to the interstellar medium iMatt eucci & Greggio|[l9860 . Be- 
cause of this, the rate of SN la explosions, and its secular evolution 
are fundamental for understanding abundances and abundance pat- 
terns in stellar systems, in galaxy clusters and in the intergalactic 
medium. Due to the high energy release of each event, the rate of 
SN la explosions also impacts on t he evolution of gas flows in El- 
liptical galaxies dCiotti et al.ll 199 lh . and on the feedback process 
which follows from an episode of star formation (SF). The rate of 
nucleosynthetic and energetic input from SN la's critically depend 
on the distribution of the delay times (DTD), the delay time being 
the time elapsed between the birth of the SN la progenitor star and 
its final explosion. In the curr ent literature, the shape of the DT D is 
still considerably debated (seeiGreggio, Ren zini & Daddil2008l and 
references therein), partly because the SN la progenitor's have not 
been unequivocally identified. 

Most notably, Type la Supernovae are crucial cosmological 
probes , e.g. having provided the first evidence for cosmic acceler- 
ation (Ries s et al]ll998l : IPerlmutter et al.lll9990 . Their use in cos- 
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mology is subject to the standardisation of their light curve, i.e. 
the ability to relate a distance independent quantity to the abso- 
lute magnitude at maximum light. For example, the Phillips relation 
links the SN la magnitude at maximum (e.g. MB. raa x) to its dimming 
in the first 15 days after maximum (Ami;, a measure of the de- 
cline rate) dPhillipd[l993l : lAltavilla et alj|2004) . While the relation 
very likely reflects different production of r adioactive 56 Ni in the 
explosion (JArnett 1982; M azzali et alj 2001), the physical mecha- 
nism which causes this diversity is not clear, nor is the reason for 
the observed dispersion of the relation, part of which could be of 
intrinsic origin. In most practical applications the standardisation 
of th e SN la light curves in volves the use of the stretch parame- 
ter fs lPerlmutter et al]| 19990 . which is related to the width of the 
light curve and the magnitude at maximum light, based on obser- 
vational data. Therefore, at present, the standardisation of the SN 
la light curves is fully empirical, and based on events at low red- 
shift: if a systematic different relation holds for high z events (either 
in the average value, or in the dispersion, or both) the cosmological 
application of SN la's as distance indicators would be called into 
question. Big, homogeneous samples of nearby events are currently 
being constructed, from which robu st correlations can be derived, 
and peculiar events single out (e.g. iHicken et al.1 120090 . but ulti- 
mately, we want to assess whether the SN la light curve stretching 
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works in the same way at high and low redshift. It is then foremost 
important to understand the diversity of SN la events, and to check 
whether this is related to parameters like age and metallicity, which 
both increase systematically from high to low redshift. 

Several observational attempts have been performed to clarify 
this issue. A relation between the 56 Ni mass and the metallicity of 
the prog enitor is expected on the basis of the nucleosynthetic pro- 
cessing dTimmes. Brown & Truranll2003l) ; this prediction has been 
studied b y trying to determ ine the metallicity of individual progen- 
itors ( e.g. Ellis et aU2008l). as well as that of the stellar population 
hosts (Gallagher et al. 2008; Howell et al. 2009). So far the results 
are not conclusive, and metallicity does not ap pear to be the most 
important parameter driving the diversity (see lHowell et al.ll2009l 
for a thorough discussion). In this respect, it is worth noting that 
a sizable fraction of mass in stars in the loc al universe be longs to 
populations that are both old and metal rich ( Renzin fcood and ref- 
erences therein); therefore also at high redshift an important frac- 
tion of SN la events come from high metallicity stars. Thus, un- 
less the metallicity dependence of the SN la light curve is very 
pronounced, the trend of its properties with redshift related to the 
change of (average) metalliticy should be weak. More significant 
are the empirical correlations found between the light curves and 
the youth of the galaxy host. Literature res ults include the aver - 
age Amis being higher in later galaxy types dAltavilla et al .1120041) ; 
brightest events occurring pr eferentially (if not only ) in hosts with 
high SF rate per unit mass (Gallagher et al. 2008); a correlation 
between the 56 Ni mass and the average age of the SN la host, 
with o nly faint events in galaxies with aver age ages older than 



3 Gyr dHowell et al1l2009l ; iNeilfe t al. 2009). The observed trends 
all point in the same direction, i.e. the younger the stellar popula- 
tion the more energetic its SN la, but the details differ. For example, 
Galla gher et al.l (2008) find that SN la's in galaxies with high spe- 
cific SF span a narrow range in V magnitude at maximum, while 
the diversit y is much more pron ounced at low specific SF. On the 
other hand, iHowell et all |2009) find a wide range of 56 Ni mass in 
events hosted by young stellar populations, which suggests t hat di- 
versity is a feature of star forming galaxies. lAltavilla et al.l d2004l) 
find a wide range of SN la light curv e properties in both early and 
late type galaxies; iNeill et alj (2009) find that faint events occur 
only in high mass galaxies, suggestive of a uniformity of SN la's in 
old stellar populations. Likely, the picture is confused because dif- 
ferent authors use different proxies to characterise the light curve of 
SN la's (stretch parameter, Am^, magnitude at maximum in some 
photometric band, ejected 56 Ni), and the age of the parent stellar 
population (light or mass averaged age, spectral indices, morpho- 
logical type). The empirical route to study this issue is further com- 
plicated by the effect of exti nction which is difficul t to correct for 
in the individual events (e.g. lNobili & Goobaj|2008l) . 

From the conject ural point of view, the claim of two popu- 
lations for SN la's I Mann ucci, Delia Valle & Pa nagia 200G), one 
(prompt) supplying events within the first 0.1 Gyr from the SF 
episode, and one (tardy) providing explosions over the whole re- 
maining age range (but with a much lower efficiency) has suggested 
the existence of two distinct channels for SN la production, possi- 
bly characterised by different light curve parameters. Since at high 
redshift the SN la production via the prompt channel must prevail, 
a systematic evolution of the Phillips relation could be present, re- 
lated to the change of mixture of prompt and tardy explosions. 

From the stellar evolution point of view, while it is generally 
agreed that SN la explosions are due to the ignition of fuel in CO 
white dwarfs (WD), the specific path to the successful event is still 
debated. The lack of hydrogen lines in most SN la spectra has tradi- 



tionally lead to exclude single stars (at the end of their evolution on 
the Asymptotic Giant branch) as viable candidates for SN la's. The 
explosion of a (naked) WD may be obtained via accretion from a 
close companion in a binary system, either as a consequence of He- 
lium detonation, when the accreted He layer reaches a critical mass 
(sub-Chandrassekhar models, hereinafter Sub-Chandra), or of Car- 
bon deflagration, when the CO WD reaches the Chandrassekha r 
mass (hereinafter Chandra) (see iHillebrandt & Niemevei] 12000), 
Different progenitor models correspond to different companions of 
the WD, which can be non degenerate (MS or post MS) or de- 
generate stars (another WD); the relative evolutionary channels are 
commonly known as Single (SD) and Double Degenerate (DD), 
respectively. Therefore, stellar evolution in a close binary system 
does provide two different channels for SN la production, and also 
two different kinds of explosion. In each case, though, the range of 
delay times can be very wide, from ~ 40 Myr (the lifetime of an 8 
M Q star, i.e. ~ the most massive CO WD progenitor) up to a Hub- 
ble time, so that each channel can accommodate prompt and tardy 
events. At the same time, a systematic variation of the exploding 
WD parameters (e.g. the range of initial mass of the primary in 
successful systems) with the delay time is expected within each 
channel, and this could account for the observed diversity, as well 
as imply correlations with the age of the parent stellar population. 

From the observational point of view, there is no clear indi- 
cation decisively supporting one channel over the oth er. The cir- 
cums tellar material detected for the event SN 2006X dPatat et all 
120071) has been interpreted as favouring the Single Degenerate 
model, while the lack of it in other events is more easily under - 
stood within the Double De generate model (see. Si mon et al .120091) . 
According to Kasen (2010) theoretical computations, the presence 
of a companion should cause an excess emission in the early light 
curve, which originates from the collision of the ejecta with the 
companion. The phenomenon, which should be detectable for ~ 
10 % of the SD events, has not been reported so far, and pro- 
vides a very interesting test of the progenitor model. Statistics 
of potential progenitors for the two models is also not conclu- 
sive. On the one hand there are many types of i nteracting binaries 
conta ining a WD plus a non degenerate star l lParthasarathv et all 
2007), like cataclismic binaries, symbiotic systems and the par- 
ticularly promising class of Supersoft X-ray sources, but it's not 
clear whether the population of good candidate SN la precursors 
among them is abundant enough to account for the observed SN la 



rate Munari & Renzinill992l ; lKenvon et alJl993l ; |Pi Stefano et al.l 
2009). On the other hand, the inventory of close double degener- 
ates from the SPY project did yield one system with total mass 
larger than the Chandr assekhar limit and clos e enough to merge 
within a Hubble time (Napiwot zki et al.ll2005l) , showing that po- 
tential SN la progenitors can be found among DDs as well. A few 
events have been reported to require a Super-Cha ndrassekhar WD 
mass dHowell et al-lhood : I Yama naka et alj|2009h : these are more 
easily explained as DD mergers. Meanwhile, single degenerate sys- 
tems containing a WD with mass clo se to the Chandrassekhar limi t 
have been discovered, like RS Oph ( Hachisu. Kato & Lunall20o"7l) 
and HD 49798/RX J648.0-4418 dMereghetti et alj2009l) . which are 
very good candidate SN la precursors. These evidences may well 
suggest that both channels are providing SNIa events with compa- 
rable efficiency. 

In this paper some general features of the DTD which are ex- 
pected from stellar evolution in close binaries are illustrated, fol- 
lowing the analytical parametrisation presented in iGreggiol d200l) 
(Paper I), and described in Sect. 2. The dependence of the DTD on 
the major parameters is considered next, focusing on the possibil- 
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ity of realising strong discontinuities which may justify the concept 
of prompt and tardy events. The systematic variation with the de- 
lay time of some properties of the progenitors (i.e. WD cooling 
times, mass of the CO WD at the start of accretion, total mass of 
the DD systems) is also examined. These two aspects impact on the 
possibility of systematic differences between low and high redshift 
samples of SN la, and are discussed in Sect. 3, separately for SD 
and DD models. Mixed (SD + DD) DTD models are presented in 
Sect. 4, where the possibility of constraining the mix from the ob- 
servations is discussed. The results are summarised and discussed 
in Sect. 5; conclusions are drawn in Sect. 6. 



2 THE ANALYTIC APPROACH TO THE DTD 

The standard evolutionary scenario to SN la explosion starts with a 
close binary system in which the primary is an intermediate mass 
star, so that the first Roche lobe overflow (RLOf) leaves a CO WD. 
Low mass stars (M < 2M Q ) do produce CO WDs, but in a close bi- 
nary they can do so only if the first mass exchange takes place after 
the He flash. Since the stellar radii at He flash are large (hundreds 
of solar radii), and since the distribution of the original separations 
Ao scales as 1/Ao, in most cases low mass primaries in close bi- 
nary systems will fill their Roche lobe (RL) while on the first Red 
Giant branch, thereby producing Helium, rather than CO, WDs. 
Therefore, reasonable limits to the primary mass (mi) in SN la pro- 
genitor systems are 2 and 8 M Q . When the secondary fills its RL 
there are two possibilities: either the CO WD accretes and burns 
the Hydrogen rich material remaining compact inside its RL, or a 
common envelope (CE) forms around the two stars. In the first case, 
the WD grows in mass and when conditions are met, a SN la ex- 
plosion occurs; in the second case, due to the action of a frictional 
drag for ce, the CE is eventuall y expelled, while the binary system 
shrinks dlben&Tutukovll 19871) . If the secondary is relatively mas- 
sive, the outcome of the CE is a CO WD + Helium star system. 
Further evolution causes the He star to expand and fill its RL, pour- 
ing Helium rich material on the companion, and the CO WD has 
again t he possibility of accreting and growing up to explosion con- 
ditions dlben & Tutukovll 19941 : IWang et alj20091) . Both these paths 
to the SN la event are SD channels, but in the first one (H-rich SD 
channel) some H is expected in the vicinity of the SN la which 
may show up in the spectrum, while in the second one (Helium star 
channel) no H-rich material should be around. 

If the CO WD does not accrete the donated Helium, a CE 
phase sets in again, leading to a close CO+CO WD system. An- 
gular momentum losses via gravitational wave radiation will even- 
tually lead the system to merge, and, if its total mass exceeds the 
Chandrassekhar limit, explosion may arise. This is the DD channel 
to SN la, as provided by intermediate mass binaries. A DD system 
may also be formed when the secondary is a low mass star, which 
fills its RL when its He core is degenerate. In this case, the CE phase 
leaves a close binary with one CO and one He WD, but, typically, 
the to tal mass of the DD system does not reach the C handrassekhar 
limit jYungelsonl [20051 ; iRuiter, Belczvnski & Frverll2009h . There- 
fore this channel provides mainly sub-Chandra events; while poten- 
tially important for the SN la's, especially at late delay times, this 
path is neglected here because sub-Chandra explosions do not pro- 
vide a good fit to light curves and spectra of the majority of SN la's 
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vide a good nt to light curves and spectra ot tne majority ot sjn la s 
i Hoefl ich et al ] | 19961 : lNugentetal.il 19971 : iHillebrandt & Niemevej 
2000). 

Upon specifications of the distribution of the initial parameters 
(m\ ,ni2, Aq), Binary Population Synthesis (BPS) codes compute the 
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Figure 1. Distribution of the delay times for SD models, Chandra (solid) 
and Sub-Chandra (dot-dashed). The DTD is normalised to 1 over the range 
0.04 to 13 Gyr. The upper axis is labelled with the value of the stellar mass 
(in Mo) evolving off the MS at the corresponding a delay time labelled on 
the lower axis. A Salpeter IMF for the primaries, a flat distribution of the 
mass ratios, and m\ > 2 (3) M Q for Chandra (Sub-Chandra) models have 
been adopted. The dashed line shows the trend of the CO WD formation 
rate from the primaries, which may represent an upper limit to the slope of 
the DTD for the MS+WD channel. This function, plotted on the logarithmic 
scale, is not normalised to 1 over the delay time range. 



evolution of a population of binaries through these complicated 
paths, determine the partition of the population in the various chan- 
nels, and the delay times of the explosions. The results are sensitive 
to the recipes used to follow the evolution, most notably the criteria 
to decide what happens when a star fills its RL. Critical prescrip- 
tions regard the occurrence of CE or not, the outcome of mass ex- 
change, the accretion modalities on the WD, whether fuel ignition 
leads to accretion induced collapse or to a thermonuclear explo- 
sion dNomoto & Iben|[l985l ; lYoon et alj|2007l ; iNomoto et ai]|2009l: 
IPiersanti et alj|2009l) . The DTD from simulations may present dis- 
continuities due to the switch on/off of one of the channels, and/or 
to an abrupt variation of the efficiency with which one channel is 
providing SN la's, as a function of the original binary parameters. 
For example, a narrow DTD, centred on a characteristic delay time, 
can be produced by suppressing all but one of the channels, and 
narrowing the range of the progenitor masses. These (putative) dis- 
continuities, and their location in delay time, follow from the pre- 
scriptions adopted in the BPS code. Many of the parameters which 
need to be specified are very poorly constrained (e.g. the efficiency 
ttcE of the CE phase which controls the degree of shrinkage of the 
system); therefore the specific realisations can be up to date, but 
remain uncertain. 

An alternative approach to derive the DTD of SN la explosions 
consists in characterising its shape using general stellar evolution 
arguments, and introducing a parametrisation of a few astrophysical 
variables believed to play an importa nt role. This approach was fol- 
lowed^ by Greggi q & Renzinil dl983h for the SD channel, improved 
in lGreggio[ 7l996) to accommodate Sub-Chandra explosions, and 
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further developed in Paper I, where analytic functions for the DTD 
are derived for both the SD and the DD channels. Compared to 
the numerical results from BPS, these functions provide a more 
flexible tool to investigate on the predictions from different pro- 
genitor models on various observables, thereby offering a way to 
constrain those astrophysical parameters which are recognised as 
crucial. On the other hand, in order to derive the analytic functions, 
some approximations and assumptions have to be made, which are 
summarised in the following. We refer the reader to Paper I for the 
detailed derivation of the analytic DTDs. 

2.1 The Single Degenerate channel 

Regardless the specific path, the delay time for SD channels is very 
close to the evolutionary lifetime of the secondary star off the MS 
(jairna)): for the H-rich SD explosion the extra time spent to reach 
the RL dimensions, and the accretion timescale during which the 
CO WD grows to ignition conditions are both much smaller than 
TrSjn.2). For the Helium star channel, some correction (of about 
10%) should be applied to account for the He burning lifetime of 
the former secondary; in practise however T n (m 2 ) fairly represents 
the total delay time also in this case. Therefore, for the SD chan- 
nel there is a one to one correspondence between the mass of the 
secondary and the delay time, so that the distribution of the delay 
times is: 

/ia(r) oc n(m 2 ) ■ \ni 2 \ (1) 

where n(m 2 ) is the distribution of the secondary masses in systems 
which provide successful explosions through this channel, and \mi\ 
is the rate of change of the stellar mass evolving off the main se- 
quence. The latter factor is a robust prediction of stellar evolution 
theory; the first factor can be derived analytically by convolving 
the distribution of primary masses with the distribution of the mass 
ratios (q = m 2 lm\) in binary systems. If all SN la's where aris- 
ing from this evolutionary path, the range spanned by m 2 should go 
from ~ 8 M Q down to < 1 M Q , to account for the wide range of 
delay times implied by SN la events happening in all galaxy types 
( Gre ggio & Renzinil 19831) . The solid line in Fig. Q]shows the DTD 
for the SD Chandra channel under this hypothesis. The two cusps 
mark the setting in of two restrictions on the parameter space lead- 
ing to explosion: the first corresponds to the minimum mass of the 
primary, which is set to 2 M Q since the exploding WD should be 
composed of CO, rather than He rich material. The later cusp is re- 
lated to the requirement of totalling a Chandrassekar mass to ensure 
explosion with a relatively low mass secondary: as the delay time 
grows long, the progressively lighter secondary ought to combine 
with a progressively heavier CO WD (hence heavier mi) to meet 
the Chandrassekhar limit. When the lower limit on m\ becomes 
larger than 2 M Q the second cusp occurs. This restriction does not 
apply to Sub-Chandra events (dot-dashed line); thus, in this case, 
the DTD presents only one cusp, motivated by the requirement of 
a minimum primary mass of 3 M Q as suggested by the models in 
IWooslev & WeaveJ < fl994h . 

A final remark concern s the WD + MS channel proposed by 
van denHeuveleIai]jl992|) : in this model the companion fills its 
RL before core H exhaustion, and the delay time could be ap- 
preciably shorter than the evolutionary lifetime of the secondary. 
The derivation of an analytic DTD is hardly possible in this case, 
since the time at which RLOf happens depends on the evolution of 
the binary separation as determined by several physic al processes , 
e.g. CE evolution, magnetic braking, tidal interactions ( llbeiill 19911) . 
However, the primary star still has to evolve and produce a CO 



WD; therefore, for a given system, the delay time will be com- 
prised between the evolutionary lifetimes of the primary and that 
of the secondary. If accretion starts soon after the formation of the 
CO WD, the distribution of the delay times will be given by Eq.Q} 
with mi replacing m 2 , i.e. the WD formation rate from the primary. 
If accretion is delayed until the secondary reaches core Hydrogen 
exhaustion the DTD will be given by Eq.lQJ. The dashed line in 
Fig. [U shows the CO WD formation rate (arbitrarily shifted) from 
primary stars with mass between 2 and 8 M Q , having assumed a 
Salpeter IMF. We may expect that the DTD for the MS+WD chan- 
nel has an intermediate slope between the dashed and the solid line, 
i.e. more skewed at the short delays compared to Eq.lfTJl, but not 
dramatically different. 

2.2 The Double Degenerate channel 

The delay time for the DD systems is the sum of the MS lifetime of 
the secondary (T n (m 2 ), hereafter t„) plus the time taken by the grav- 
itational wave radiation to get the DD system into contact (t G w)- 
The latter time delay depends on the separation (A) and masses of 
the components of the DD system at birth; if both ;«[ and m 2 are 
intermediate mass stars (> 2M G ) this delay can be approximated 
as tgw = 0.6 A 4 M^p Gyr, where Mdd is the total mass of the DD 
system, and both A and M DD are in solar units (Paper I). For what 
explained above, this limit on the masses of both binary compo- 
nents corresponds to considering CO + CO DD systems; besides, 
the limit m 2 > 2 M Q implies a maximum delay T n of 1 Gyr. The de- 
lay tgw, plotted in the left panel of Fig. [2] is very sensitive to both 
A and Mdd, and in principle it ranges from virtually zero, for the 
most massive / closest systems, up to infinity for the least massive 
/ widest ones. Fixed limits apply to Mdd, from 1.4M Q (Chandra 
explosions) to ~ 2.4M Q , when both components in the primordial 
system are of 8 mJJ More elusive is the allowed range for A, al- 
though it must be smaller than ~ 4R in order to provide events 
within a Hubble time. 

Typically, the separation of the DD systems at birth is the re- 
sult of evolution through at least 2 RLOf episodes. Given the diffi- 
culty of describing this complex hydrodynamical phase, the results 
of the application of the various recipes in the literature are un- 
certain. In this respect, it is important to notice that a substantial 
shrinkage is necessary to account for SN la events within a Hubble 
time from DD systems, since the first RLOf occurs when the sepa- 
ration A ranges from several 10 R B (to avo id premature merging) , 
to several 100 R Q (to get RLOf at all, see llben & Tutukov|[l9871) . 
Therefore, the prescription used in BPS codes to determine the de- 
gree of shrinkage should be accurate at a level of ~ 1 % to yield 
correct results on an individual binary. 

The sensitivity of the ratio A/A to the commonly adopted lit- 
erature descriptions of the CE phase is explored in Paper I. It ap- 
pears that when the ctce formalism jWebbink|[l984l) is applied to 
two mass exchange phases, the ratio is very small, so that only ini- 
tially wide systems manage to populate the relevant range of sep- 
arations (up to a few R ). Conversely, if the fi rst of the two CE 
phases is regulated by the lNelemans et all J20CXE) (angular momen- 
tum loss) recipe, the ratio A/A covers a ~ 2 orders of magnitude 
range, depending on the masses of the components. Moreover, in 
the former case, the heavier the original binary the lower the A/A 

' The remnant of an 8 M e star is ass umed to be a CO WD of 1 . 2 M , 
consis tent with the empir ical relation in Williams, Bolte & Koester j2004l) 
and in lKalirai et alj2008l 
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Figure 2. Left: the delay time due to the gravitational wave radiation, as given by the labelled (approximate) relation, as function of the initial separation of 
the DD systems, and for Mdd ranging from 1.4 to 2.4 M Q in steps of 0.1 M Q . Right: the distribution of the delay times for DD models, CLOSE (solid) and 
WIDE (dot-dashed), in the same units as in Fig. Adopted parameters include a Salpeter IMF for the primaries, a flat distribution of the mass ratios, t„ < 1 
Gyr, and flat distribution of the separations A. The DTDs for the SD models are plotted in red for comparison. 



ratio: this is because in the ci C e formalism the shrinkage at each CE 
phase is proportional to the binding energy of the envelope of the 
donor, which is larger when the donor is more massive. Thus, the 
larger m 2 , the shorter r„ and the shorter t gv/ : this skews the DTD 
towards the short delay times. When the Nelemans' formalism is 
applied, instead, no strong relation appears to emerge between the 
initial masses of the binary components and the shrinkage, so that 
t„ and r G w are likely to be uncorrelated, as well as M DD and A. 

On this premises, analytic distributions of the delay times for 
the DD channel have been derived in Paper I in two flavours: the 
DD-CLOSE model, in which at each m 2 the delay time varies be- 
tween a minimum (=T n (m 2 )) and a maximum delay, whose value 
increases as m 2 decreases. In this scenario, a wide portion of the ini- 
tial parameter space (m 1 ,m 2 ,A ) maps into short delay times, while 
the long delay times are populated with only the less massive bina- 
ries. The DD-WIDE model, instead, assumes that A and Mdd are 
independent variables, and their distributions concur to shape the 
distribution of the gravitational wave radiation delay. The variables 
t„ and tqw are also assumed to be basically independent, except 
for a correlation between m 2 (i.e. r n ) and M DD , to account for the 
fact that the heavier the original secondary, the more massive the 
DD remnant. As a result, in the DD-WIDE channel the upper limit 
to t gw is greater then the Hubble time irrespective of the mass of 
the original binary, and the population at long delay times of the 
corresponding DTD is enhanced. Simple power laws are assumed 
to describe the distribution of r GW (CLOSE DDs) and A (WIDE 
DDs), with parametrised exponents respectively called /?„ and j8 a . 
The shape of the distribution of Mdd is found to be relatively unim- 
portant. 

Admittedly, the above approximations and the assumptions 
are rough and do not punctually describe the outcome of the evo- 
lution of close binaries, SN la progenitors. However, the range 



spanned by Mdd and A for events within a Hubble time are small, so 
that the assumption of smooth distributions of these parameters is 
likely appropriate. The main features of the curves are determined 
by the gravitational wave radiation clock, which implies that the 
distribution of t G w is skewed at the short end, unless all successful 
DDs are bom very wide. This characteristic can be judged by pro- 
jecting a flat distribution of A into the corresponding distribution 
on t gw in Fig. [2] left panel. 

The right panel of Fig. |2] shows the DTD for CLOSE DDs 
and WIDE DDs, normalised to 1 over the range of delay times up 
to 13 Gyr. The distributions can be viewed as a modification of the 
SD case, where systems with given t„ populate a range of delay 
times longer than this due to the extra delay r GW . The WIDE DD 
case is the least populated at the short delay times, the CLOSE DD 
the most, since in this case the great majority of systems have a 
short t G w- Notice that delay times in excess of 1 Gyr are populated 
only with systems whose t G w is longer than some limit, since by 
construction a maximum r n of 1 Gyr is imposed. This is the rea- 
son for the presence of the cusp in the DTDs: at this delay time 
a fixed limit on the allowed parameter (r n ) space is hit. Although 
this feature in the analytic DTDs is purely numerical, we do ex- 
pect an abrupt change in the distribution at a delay time equal to 
the MS lifetime of the least massive secondary in successful SN la 
progenitors. 

Inspection of Fig. [2] shows that the various models provide 
similar curves: the DTD features an early maximum and most 
events occur within ~ 1 Gyr from formation. At later epochs the 
DTD exhibits a steady decline with a pronounced downturn for 
the SD Chandra model. All distributions accommodate prompt and 
tardy events, but with a substantial continuity over the whole range 
of delay times; this follows from the assumption of continuous dis- 
tributions over the relevant parameter space of masses and sepa- 
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rations of SN la progenitors. In the next section these curves are 
compared with theoretical results in the literature. 

2.3 Analytic DTDs and BPS results 

In deriving the analytic DTDs described above one potentially im- 
portant variable is neglected, i.e. the separation of the primordial 
binary Aq. In order to provide a successful event, Ao ought to range 
between a minimum and a maximum value, respectively to avoid 
premature merging, and secure interaction at all. As long as these 
limits do not depend on the component's masses, the shape of the 
DTD is not affected, and the limitation just reflects into the over- 
all likelihood of the SN la event from a stellar population. How- 
ever, the minimum and maximum initial separations tend to in- 
crease with the mass of the primary, because more massive stars 
have larger radii at all key points along evolution. Figures 1 to 4 
in Iben and Tutukov (1987) well illustrate this point, also showing 
that, in general, the range of separations for a successful interac- 
tion widens as the primary mass grows. It the distribution of A 
were flat, the more massive binaries would have then greater chance 
to provide SN la's; on th e other hand, closer systems should be 
more abundant jAbTl983). so that the more massive binaries also 
have a higher likelihood of premature merging. Therefore, it is not 
clear how the efficiency of SN la production from a binary of given 
(ra l ,ra 2 ) varies with the system separation. The comparison of the 
analytic DTDs with results from BPS helps understanding whether 
or not the analytic functions capture the most important astrophys- 
ical parameters. 

In Paper I the ana l ytic D TDs are compared to the BPS results 
of l^ngelson & Livio (2000): consistency is found for the Sub- 
Chandra and for the DD-Chandra cases, when accounting for the 
limits on mi and mi adopted by the authors. The comparison be- 
tween the DTDs of the SD-Chandra channel is hampered by the 
narrowness of the range of secondary masses in systems which pro- 
vide S N la through the H-rich SD path in the I Yungelson & Livid 
(2000) simulations. Actually, the conditions for stable mass trans- 
fer and/or a successful growth to the Chandrassekhar limit of the 
H-accreting CO WD are often found to hold only within a very 
restricted range of m 2 . The analytic DTDs do not incorporate con- 
straints on the mass transfer rate, but the mass range covered by the 
components of successful systems is to be considered as a parame- 
ter. 

lHachisu, Kato & Nomotcl ( 120081) present the DTD for the SD- 
Chandra scenario taking into account the effect of a radiative 
wind from the accreting WD w hich stabilises the mass transfer 
( lHachisu. Kato & NomotdT l996). favouring the WD mass growth. 
Although their computation is not a BPS simulation, the DTD is 
obtained by integrating the distributions over the appropriate range 
in the parameter space (mi, m 2 , A ) which leads to successful ex- 
plosions, according to their modelling of the evolution of CO WDs 
accreting from a MS and from a RG companion. It is instructive 
to compare their derived DTD with the Paper I anal ytical distri- 
bution s. Let's first consider the WD+RG channel: lHachisu et ajj 
(2008) find it active for m 2 < 3 M G , and indeed the correspond- 
ing DTD is populated at r> 0.4 Gyr; its DTD is marked by two 
cus ps, at 1.6 and 12.6 G yr. Given the wide bins of delay time in 
the lHachisu et alj j2008l) computation, the position of the cusps are 
consistent with those on the solid line in Fig. [T] For the same choice 
of IMF (Salpeter) and distribution of mass ratios (flat) as in Fig. [T] 
the slope of the DTD between the two cusps (in logarithmic units) 
is -1.4, perfectly consistent with the analytic DTD. The WD+MS 
channel is found active for 2 < m.2/M Q < 6 and, in fact, the relative 



DTD goes to zero in the vicinity of r= 1 Gyr; the slope is — 1.4, 
much steeper than the slope of the WD formation rate (dashed line 
in Fig. [TJ in the same range of delay times, which is — 0.5. This 
may signal that the distribution of the separations is important in 
shaping the DTD of the WD+MS channel, favouring SN la events 
in the more massive binaries. 

The DTD from the WD+MS ch annel has also been pre- 
sented bv lHan & Podsi adlowski (2004), as derived with their BPS 
code. The resulting DTD is a rather narrow, almost flat distri- 
bution, within limits which depend on the parameter for CE ef- 
ficiency. A similar shape is found for the WD+MS channel in 
iMeng. Chen & Hani d2009l) . The limits are consistent with the MS 
lifetime of the most and least massive secondary in SN la pro- 
genitors through this channel, which are found confined between 
~ 3 and ~ 2 M e . T hese DTDs are quite different from those in 
lHachisu et al.l d2008h . possibly because of the use of a different 
technique and of different details in the assumptions. Therefore, the 
shape of the DTD for the WD+MS channel seems rather sensitive 
to the adopted approach. 

DTDs based on BPS c al culatio ns have been presented 
by iBelczvnski. Bulik & RuiteJ ( 120051) and more recently by 
Ruite r et al.1 ( 120091) . Since the former models include Sub-Chandra 
explosions, which actually provided the majority of events 
(Ruiter ct ail 120091) . their comparison with the analytic DTDs is 
not practical. More simi lar are the SN la p roduction channels in 
the models presented in iRuiter et"al1 d2009h . which include only 
Chandra exploders, from DD progenitors (DDS), H-rich SD Chan- 
dra progenitors (SDS), and CO WDs accreting Helium via stable 
RLOf, either from a non degenerate or from a degenerate compan- 
ion (AM CVn). The DTD for the DDS channel is very similar to the 
analytic function, with a logarithmic slope of =-1. Notice that when 
the efficiency of the CE phase is reduced (acE=0.5) the DTD of the 
DDS channel is found steeper, as a consequence of the smaller or- 
bital separations of the DD systems when they emerge from the 
CE. The same effect can be obtained by varying the parameter f} g 
or /3 a in the analytic DTDs. The comparison with the BPS results 
for the SD channel is more comp licated, due to its sp litting into the 
two classes SDS and AM CVn in lRuiter et"aT1 d2009l) . However, the 
BPS DTD seems much flatter than the analytic DTD in the range of 
delay times between 3 and 6 Gyr, and lacks the dramatic drop at late 
epochs. Since the clock of the explosion is the evolutionary lifetime 
of the secondary (also in the BPS simulations), it seems that the en- 
velope mass of the secondary has little influence on the efficiency 
of the production of Chandra expl osions in the SDS channel, in the 
simulations by Ruiter et al. (20(3). 

To summarise, the comparison of the analytic curves with 
other results in the literature shows that the limits on the compo- 
nent's masses play a crucial role in shaping the DTD, while the 
effects of the limits on A are much less evident, except perhaps 
for the WD+MS progenitors. While the analytic DTDs for the DD 
models match well the BPS results, the DTD of SD models turns 
out appreciably different. It seems however that diverse realisations 
are obtained by different authors, both for the range of m 2 of SN la 
progenitors, and for the DTD slope. These differences likely reflect 
the sensitivity of the results of the systems' evolution to details in 
the prescriptions adopted for the simulations. 



3 EFFECT OF PARAMETERS 

By construction, the analytic DTDs have a built in parametrisation 
of some astrophysical variables related to the binary population and 
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Figure 3. Effect of the distribution of the mass ratios on the mass distribution of the secondaries (left) and on the DTD (right) of SN la from the SD-Chandra 
channel. All curves assume (m\ jmi) *S 8M Q , an d a Salpeter distribution for the primary masses. Red lines assume f(q) oc q? with y = (solid) and 1 (dashed). 
Black lines assume Eq.(2} (solid) and (3) (dashed). The cumulative distributions of q are shown in the inset (left panel). The dash-dotted blue line in the left 
panel shows a Salpeter distribution: with respect to this line, 11(1112) is underpopulated at both the high and the low mass end, because of the requirements on 
the binary system SNIa progenitor. The deficiency at the high mass end is due to the upper limit of 8 M Q on the primary mass; at the low mass end, instead, it 
is due to the lower limit on the mass of the WD. 




Figure 4. Distribution of the secondaries (left) and on the DTD (right) of SN la from the SD-Chandra channel for a flat f(q). All curves assume (mi ,mi) 
< 8M Q , but differ in the value of other parameters: the solid and dotted lines assume a minimum m\ of 2 and 5 M B respectively, no lower limit on 1112, and 
e=l; the dashed line assumes the same limits as the solid line, but e = 0.5; the dot dashed lines assumes a minimum ni2 of 4.5 M Q . The line encoding is the 
same in the two panels. The top axis in the left panel is labelled with the MS lifetime of the evolving secondary, in Gyr. 
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to the progenitor's model. In addition, they are constructed by re- 
lating the delay time to the exploding system, so that systematic 
variations of the properties of SN la progenitors with the delay time 
can be envisaged. In this section these issues are illustrated for the 
Single and the Double Degenerates. 

3.1 Single Degenerates 

The curves in Fig.Q]are constructed upon specification of: the slope 
of the mass distribution of the primaries (a), the distribution of the 
mass ratios (f(q)), the mass range of the primaries and of the sec- 
ondaries in successful systems. An additional parameter is the ac- 
cretion efficiency e, namely the fraction of the envelope mass of 
the secondary which is actually accreted by the primary. The solid 
curve in Fig.Q]assumes a = 2.35 (as all other DTDs in this paper), a 
flat distribution of the mass ratios, the maximum possible range for 
m t and m 2 , and that the whole envelope of the secondary is avail- 
able for the WD to grow up to the Chandrassekhar limit (e = 1). 
The same curve is reproduced as a solid line in Fig. [4] The de- 
pendence on a and f(q) have been illustrated in Paper I: basically, 
any combination which favours low mass secondaries leads to a 
larger fraction of late epoch events. In Paper I it is also shown that 
when using the alternative algorithm adopted in Gregg io~& Renzinil 
Jl983h to derive n{m{), the fraction of low mass secondaries is en- 
hanced, for a given IMF and distribution of the mass ratios^. 

Empirical deteminations of the distribution of the mass ratios 
in binaries are subject to con siderable observational biases (see 
e.g. iKouwenhoven et alj|2007h . Most binary population synthesis 
computations assume a flat distribution; on the other hand, there is 
growing evidence for a prominent population of twins among close 
binaries, w hich could enhance the fractio n of systems with short 
delay times dPinsonneault & Stan ek (2006) and references therein). 
Fig. H] shows the distribution of the secondary masses and the en- 
suing DTDs under various assumptions for fiq). In particular, the 
two black lines assume a distribution of the mass ratios made of 
two components: the solid line adopts 

f(q) = 0.63 + 2.42 x (g/0.95) 10 (2) 

which is a flat distribution up to q - 0.8, plus a prominent 
population of twins. This fu nction roughly fits the distribution in 
IPinsonneault & Stanekl ^20061) for q > 0.5. The dashed black line 
adopts 

f{q) = 6.9q 2 - 6.9^ + 2.15 (3) 

which roughly reproduces Trimble (2008)'s results, with a strong 
component at the lowest and highest values of the mass ratio, and a 
minimum at q = 0.5. Clearly, the distributions with a high propor- 
tion of twins (the two dashed lines and the solid black line) provide 
DTDs which are relatively more populated at the short delay times. 
The effect however does not appear dramatic, nor does a two com- 
ponents f(q), as in Eqs. (O and l[3}, reflect into a two components 
DTD. All following computations adopt a flat distribution of the 
mass ratios, an option which maximises the fraction of late epoch 
explosions. This will be taken into account when discussing the re- 
sults. 

Similar to Fig. [3] Fig. [4] illustrates the effect of varying the 

2 Notice that the Greggio & Renzini 1 1983) formalism assumes power law 
distributions for the total binary mass Mb and for the ratio of ni2 to Mb , so 
that the same IMF slope and distribution of the mass ratios have different 
meanings in the two approaches. 



other parameters. The dashed curves assume the same limits for 
the masses of the components as the solid curve, but e = 0.5 for 
all systems. The possibility of accreting only half of the envelope 
of the secondary has a dramatic impact on the number of systems 
with a low mass secondary which manage a Chandra explosion; 
correspondingly, the DTD drops rapidly off at delay times slightly 
over a Gyr. This makes it difficult for the SD channel to account for 
a sustained SN la rate at long delay times: for example, at an age of 
1 3 Gyr, the stars evolving off the MS for a single burst stellar pop- 
ulation have a mass of ~ 0.9M o . With a core mass of 0.3 M , only 
0.6 M are available for accretion onto the companion WD, and, if 
6 = 0.5, the Chandrassekar limit can be met only starting with a CO 
WD heavier than 1.1 M . Therefore, the range of suitable primary 
masses is exceedingly small (7 < m\ < 8). The dotted curves in 
Fig. |4] show the effect of implementing a more severe lower limit 
to m\ : 5 M instead of 2 M , corresponding to SN la's arising only 
from CO WDs more massive than ~ 0.9 M e . In this case, 11(1112) 
gets substantially deprived of systems as soon as m 2 drops below 5 
M Q , and the DTD gets depressed at the corresponding MS lifetime 
(- 0.1 Gyr); however the distribution remains wide, without strong 
discontinuities. A strong discontinuity may instead be obtained by 
implementing a dramatic change in SN la production at some value 
of ni2, due to the one to one correspondence between delay time 
and mass of the secondary in the SD model. In particular, a lower 
limit on m 2 has the effect of vanishing the DTD at delays longer 
than the MS lifetime of such stellar mass.The dashed-dotted lines 
in Fig. [4] show the effect of adopting a lower limit of 4.5 M G to 
the secondaries in SN la progenitors, mimicking what happens for 
the SD Helium star cha nnel which require s massive secondaries. 
Notice that, according to lWang et al J d2009h . the secondary mass in 
such SN la progenitors ought to be larger than 5.6 M , implying a 
drop off in the DTD at ~ 80 Myr. 

I turn now to consider some properties of the exploding WD 
which are predicted to vary sys tematically wi t h the delay time 
in the context of the SD model. iLesaffre et alj |2006) presents a 
systematic study of the sensitivity of ignition conditions for H- 
rich Chandra SD exploders on various properties of the progenitor. 
These authors find that the more massive and/or the cooler the CO 
WD when accretion starts, the higher the central density at ignition. 
While it is not clear whether this leads to a larger or lower mass of 
synthesised 56 Ni, a systematic variation of ignition density could 
translate into a systematic variation of the properties of the SN la 
explosions. It is thus interesting to analyse how the initial mass of 
the CO WD and its cooling time vary with the delay time. 

The top panel in Fig. [5] shows the range for the initial WD 
mass of SN la progenitors as function of the delay time, for a few 
choices of the relevant parameters. The limits on the WD mass fol- 
low from the limits on m\ through the adopted initial-final mass 
relation (M wd = 0.6 + 0.1 x (M - 2.) as in Paper I). The maximum 
value of the primary mass is 8 M Q irrespective of the delay time; 
the lower limit is given by mj (function of t), or higher, either be- 
cause a fixed limit on m l sets in (2 M Q for the solid line, 5 M e 
for the dotted line), or, at late delay times, because of the need of 
securing a Chandrasekhar mass at explosion. Correspondingly, the 
range of initial mass for the CO WD is small at short delays (only 
massive binaries have had time to evolve), but it rapidly increases 
with r; at intermediate delay times the range of initial WD mass is 
maximum, but at late delays it shrinks, because a massive primary 
is required to accomplish a Chandra explosion. The latter regime 
sets in earlier for a less efficient accretion (e.g. the dashed curve). 

The bottom panel of Fig. [5] shows the range of cooling times 
of the CO WD, i.e. the time elapsed between its the formation and 
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the start of the accretion phase. In this time interval the WD cools, 
becomes progressively more degenerate, a nd other processes like 
C and O separation and crystallisation (e.g. Sala ris et al . 2000) set 
in. The path towards explosion, and the explosion outcome criti- 
cally depend on the initial structure of the CO WD, or, in turn, on 
the cooling time. For example, the C and O relative diffusion may 
affect the amount of synthesised 56 Ni, hence the light curve. On av- 
erage, the cooling time increases with the delay time, and its range 
reaches a maximum at t- 1 Gyr, to become smaller at longer r. 
Most noticeably, the cooling time is very long, up to several Gyr, 
for most values of delay time. 

To summarise, some dispersion in both parameters is present 
at all delay times, but for early events this dispersion is more pro- 
nounced. Qualitatively, then, we may expect more diversity in SN 
la events from younger stellar populations. On the other hand, at 
short delays the CO WD at the beginning of the accretion phase 
is (on average) more massive and younger; thus, due to compen- 
sating effects, the trend of ignition density with delay time could 
be relatively weak. Conversely, at long delay times the mass and 
cooling age of the CO WD both increase with r, so that the ignition 
density should increase systematically. It is important to notice that 
very long cooling time s, of few to several G yrs, characterise most 
of the events, whereas iLesaffre et all (2006) models include cool- 
ing ages of ^ 1 Gyr at most. As the delay time increases, the CO 
WD has time to become more and more degenerate before the start 
of the accretion phase. The evolution of the WD structure during 
the accretion phase, as well as the outcome of C ignition under ex- 
tremely degenerate conditions are important questions which need 
more attention than currently present in the literature. 



3.2 Double Degenerates 

Due to the relatively restricted mass range of the binary compo- 
nents (i.e. from 2 to 8 M ) in the DD model, the DTD is not very 
sensitive to the IMF slope and to the distribution of the mass ratios, 
whereas the main parameters controlling its shape are the maxi- 
mum delay due to the MS evolution of the secondary (t„_ x ), and 
the shape of the distribution of the separations. The parameter T n x 
is the evolutionary lifetime of the least massive secondary in suc- 
cessful SN la progenitors: the heavier this limit is, the earlier the 
position of the cusp in the DTD for both CLOSE and WIDE DDs. 
The shape of the distribution of the separations of the DD systems 
controls instead the overall slope of the decline as the delay time 
grows: if most systems are born close, most events have a short r GW 
and the DTD declines faster. Fig. |6]shows these dependences; no- 
tice that the DTDs are not normalised to 1 for illustrative purposes. 
The three values of r nx which characterise the curves in the left 
panels correspond to a minimum mass of the secondary in success- 
ful systems of 2 (solid), 2.5 (dashed) and 3 M Q (dot-dashed): for 
normalised DTDs, the fraction of early events is higher when this 
minimum mass is larger. The curves in the right panels all adopt 
a minimum m-x of 3 M G , but different distributions of the separa- 
tions, as described by the parameters /? a and /?„ for the WIDE and 
CLOSE cases respectively. Although the two parameters are de- 
fined in a different way (see Sect. 2.2), the smaller their values the 
larger the fraction of systems with small separation A, and the three 
curves correspond to a flat distribution (solid), a distribution more 
populated at small (dot-dashed) and at large (dashed) separations. 
For all the values of the /? parameters considered here, the DTD de- 
clines towards long delay times, but the slope is milder when f3 a or 
f}„ are larger, since in this case more systems have a long gravita- 




log DT (yr) 

Figure 5. SD model: systematic variation with the delay time of the mass 
range (top) and of the cooling time range (bottom) of the CO WD when ac- 
cretion starts. The different curves refer to different choices of the minimum 
m\ and of the parameter e with the same line encoding as in Fig. [4] 

tional wave radiation delay. Consequently, the higher /? a or /?„, the 
larger the fraction of events at late epochs. 

Some other parameters have been fixed to derive the DTDs 
shown in Fig. [6] the minimum evolutionary lifetime of the sec- 
ondary (r n o) and the minimum t G w ( t gw,o) m successful systems. 
The adopted values are (T n o,T GWO ) = (0.04,0.001) Gyr, and cor- 
respond, respectively, to assuming that (i) the most massive sec- 
ondary in SN la progenitors is an 8 M G star, and (ii) that the mini- 
mum delay due to gravitational wave radiation is negligibly smallU 
On the other hand, the upper limit to mi could be lower than 8 M Q , 
if, e.g., systems with massive secondaries feed the SD, rather than 
the DD channel. Likewise, a significant lower limit to t G w could 
be realised in nature, if, e.g., DD systems which are born too close, 
or are too massive, undergo an accretion induced collapse rather 
than a SN la explosion. Both alternatives deprive the DTD of early 
events, while maintaining the population at long delay times. Fig. [7] 
shows the effect of varying separately t„ and r GW for the CLOSE 

3 The value of 1 Myr is imposed for numerical reasons. 
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Figure 6. Effect of the main parameters on the DTD of DD CLOSE (upper 
panels) and DD WIDE (lower panels) models. The left panels show the 
dependence on the maximum evolutionary lifetime of the secondary: r n x = 
1 (solid), 0.6 (dashed) and 0.4 (dash-dotted) Gyr, for the labelled values of 
/S a and yS„ . The right panels show the dependence on the distribution of the 
separations as labelled. Notice that the DTDs are not normalised to unity 
for better illustration. 



and the WIDE DD models: some DTDs assume a nominal tgw.o 
coupled to t„ o=0. 1 and 0.15 Gyr (corresponding to a minimum m 2 
of 5.2 and 4.3 M B respectively); some DTDs assume the default 
t„_o=0.04 Gyr coupled to tgw,o=0.1 and 0.2 Gyr. The latter limits 
correspond to systems with a separation of 1 R Q and masses Mdd=2 
and 1 .4 M Q respectively. In all cases the DTD is deprived of events 
with delay times shorter than (r n ,o + tqw.o), but the overall shape is 
not much affected, though for the WIDE DD case a higher r G w,o 
spreads the DTD toward longer delay times. It is also worth noting 
that the cusp in all DTDs is actually located at t=t„. x + tqw.o: in- 
deed, when the minimum gravitational wave radiation delay is non 
negligible, the limit t„ < r„ x corresponds to a discontinuity in the 
parameter space for r located at T n-X +T GW0 . 

For the DD channel, at each delay time the events come from 
systems with a range of t„ (e.g. T n ,o < t„ < r - tqw.o) and a cor- 
responding range of t G w - Therefore, systems with different masses 
of the components and different separations end up exploding at 
the same delay time, so that the occurrence of a significant discon- 
tinuity in the DTD at a particular delay time is unlikely. Similarly, 
there is no one to one correspondence between the delay time and 
the mass of any of the WD components; however there is a trend 
of high mass binaries providing preferentially explosions at early 
delays. 

The diversity of the SN la light curves could be due to the 
range of total mass Mod if, e.g., the amount of 56 Ni produced in 
the event was correlated to the total mass of the exploding system. 
For example, the extra mass (M DD - 1.4) could act as a tamper 
in the explosion, i.e. inertially containing it a little longer, and lead 
to a higher production of 56 Ni. Although this possibility has not 
been explored with theoretical models, it is worth noting that the re- 
markable event SNLS-03D3bb was both over-luminous and Super- 
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Figure 7. Effect of varying the minimum evolutionary lifetime (r„ o) and 
gravitational wave radiation delay (tqw.o) in SN la progenitor systems 
for DD CLOSE (top) and WIDE (bottom) models. The line-types encodes 
the values of (r„,o,T G w,o) = (0.04,0.001) (solid), (0.1,0.001) (dot-dashed, 
red), (0.15,0.001) (dotted,red), (0.04,0.1) (short dashed), (0.04,0.2) (long 
dashed). All curves assume a minimum mi of 2.5 M Q , and a flat distribu- 
tion of A. 

Chandrassekhar (Ho well et alj2006h . Fig.[8]illustrates the expected 
trend of the range of M DD with the delay time. The lower limit rep- 
resents the less massive systems which have had time to evolve at 
a given r, i.e. binaries in which m 2 is the evolutionary mass with 
lifetime equal to r, and m\=ni2. For a negligible value of tgw.o, 
some of these systems merge at a delay time t. If the close binary 
evolution ends up with a wide range of separations of the DD sys- 
tems irrespective of the components' masses, the maximum Mdd 
is the remnant of the (8+8) M Q systems, i.e. M DD = (1.2+1.2) M G , 
independent of t. It will be in fact always possible to find a massive 
system wide enough, so that the gravitational wave radiation delay 
equals T-r n . This upper limit could be appropriate for the WIDE 
DD case (dashed line). However, massive binaries, which end up 
with massive remnants, also undergo a severe shrinkage, especially 
in the CLOSE DD case. Therefore, high M DD will be associated to 
short r G w, and short total delays in general; this introduces a trend 
for the upper limit to Mdd with r. The solid line in Fig.[8]illustrates 
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Figure 8. Illustration of the dependence on delay time of the range of DD 
systems total masses for the WIDE (dashed) and the CLOSE channels. For 
the latter, the upper limit shown as a solid (dash-dotted) line refers to a 
minimum secondary mass of 2 (3) M Q in SN la progenitors. A maximum ni2 
of 8 M Q and a negligible lower limit for tqw have been adopted. For delay 
times longer than ~ 0.35 Gyr the lower limit is given by the Chandrassekhar 
mass. 



this effect, and has been computed by applying the a C E recipe to 
two consecutive RLOf episodes. The exact location of this upper 
limit depends on various parameters (e.g. «ce, maximum separa- 
tion of the interacting systems, initial-final mass relation), but the 
general trend is realistic. In fact, at each delay time, the maximum 
M D d is realised in systems with maximum remnants from m! and 
ni2 which merge within tqw = t — T n (ni2). Since tqw is shorter in 
more massive systems, due to its explicit dependence on M DD , as 
well as to the higher shrinkage, the above condition is met for cer- 
tain values of m\ and mi below their maximum values of 8 M Q . The 
calculation shows that, as the total delay time grows, the condition 
above is first met with the maximum m\ (8 M Q ) in combination 
with a progressively less massive secondary, until the latter reaches 
the lower limit on m 2 (e.g. 2 M G ) adopted for the DD progenitors. 
At delay times longer than this the maximum Mdd is realised in 
systems with m 2 = 2 M Q and a progressively lighter maximum pri- 
mary. The transition between the two regimes is marked by the cusp 
in Fig! 

To summarise, for the WIDE DD case we expect a wide range 
of M DD at any delay time, and virtually no trend for the average 
mass of the exploding WD with the age of the parent galaxy. In- 
stead in the CLOSE DD case we expect a definite trend of decreas- 
ing average M DD with increasing delay time, i.e. with increasing 
average age of the parent population. In addition, the range of Mdd 
in SN la progenitors should become narrower in old systems, where 
the events should basically come from M DD close to the Chan- 
drassekhar mass. 



4 MIXED MODELS 

As shown in the previous sections, the different channels for SN 
la production are in principle capable to provide both prompt and 
tardy events, as required by the data. Moreover, for all channels 
the DTD favours early events, so that the SN la rate per unit mass 
in young stellar systems is expected higher than in old ones, again 
as required by the observations. Although the fraction of early ex- 
plosions is different for the different channels, and within the same 
channel it depends on the astrophysical parameters which charac- 
terise the progenitor systems, each model is compatible with the de- 
creasing rate per unit mass going from late to early type galaxies. 
On the other hand, there's no real need to view the various mod- 
els as mutually exclusive, and both SD and DD paths could well 
be realised in nature. From the stellar evolution point of view, the 
SD channel has more difficulty to produce Chandrassekhar explo- 
sions at very late epochs past star formation, because of the small 
mass reservoir from the donor; the DD channel instead, can pro- 
vide events at delay times exceeding the Hubble time due to the 
possibility of realising very long tqw- Therefore, rather than con- 
sidering the two channels as mutually exclusive, it seems wiser to 
construct a mixed model, in which both SD and DD events concur 
to the SN la explosion. Among the many possible combinations, 
I consider two extreme hypothesis: the Solomonic mix, in which 
both SD and DD channels contribute half of the total events within 
13 Gyr; the Segregated mix, in which the SD explosions provide 
the early events, while the DD channel feeds all the remaining de- 
lay times. Although many other combinations are possible, it is in- 
structive to see how these choices, which maximise the contribution 
of the two channels, impact on the DTD, and explore how the mix 
may be constrained by the observations. 

Fig. [9] shows the Solomonic (left) and the Segregated (right) 
mixtures obtained with the particular choice of parameters for the 
individual channels as described in the caption. The Solomonic mix 
is characterised by a DTD which is very similar to each of the SD 
and DD models used to construct it; the only significant feature 
is the intermediate slope of the decline at very late delays, where 
the number of events is kept relatively high due to the contribu- 
tion from the DD channel. This characteristic is relevant for the 
SN la rate in the oldest galaxies, as will be illustrated later on. At 
each delay time, the ratio between DD and SD explosions varies, 
but it keeps within a factor of 2 over almost the whole range of t: 
only for r > 9 Gyr does the DD channel definitely prevail. This 
kind o f mixture and its evo lution over cosmic time has been stud- 
ied in lGreggio et al.1 J2008). In general, for the Solomonic mixture, 
the majority of SN la events in early type galaxies should come 
from DD systems with M DD close to the Chandrassekhar mass. 
The events in late type galaxies, instead, should come from both 
channels and may include super-Chandra explosions from DD pro- 
genitors. 

The Segregated mix assumes that systems with m 2 > 5M Q 
feed the SD channel, while systems hosting less massive secon- 
daries undergo CE evolution at the second RLOf, to become DDs, 
hence the discont inuity at 0.1 Gyr. This combination mimics the 
DTD proposed bv lMannucci et al.l d2006h to explain the high SN la 
rate in radio loud Ellipticals as due to a recent SF episode in com- 
bination with a high fraction of prompt explosions. Actually, if this 
were the case, the excess of 7 la events in the radio loud Ellipticals 
would imply that ~ 14 CC SNe should have been det ected in the 
same sample, while none was observed (Greggio et al. 200§|). It is 
anyway instructive to consider this DTD as an example of a distri- 
bution with a high fraction of events at very short delay times; here 



12 



Laura Greggio 




Figure 9. DTDs for mixed models. Left: the DTD of SD Chandra explosions in Fig.[T](dotted line) is combined with the DTD of a DD-CLOSE model with 
minimum secondary mass of 2.5 M Q and a flat distribution of A (dashed) to construct the Solomonic mixture, in which 50 % of the total events within 13 Gyr 
come from either channel. Right: the events at delay times shorter than 0. 1 Gyr come from SD explosions, while those at longer delay times come from the 
DD-CLOSE (solid) or DD-WIDE (dot-dashed) evolutionary channels, both with minimum im of 2.5 M Q and a flat distribution of A. The events from the SD 
channel amount to 30 % of the total within 1 3 Gyr. 



a fraction of 30 % prompt events was chos en following the indica- 
tions from chemical evolution models by Matt eucci et al.l |2006). 
In this option, all events in early type galaxies come from the DD 
channel, while in galaxies with ongoing star formation we get a 
mixture of SD and DD explosions. In most cases, the two channels 
will be providing events with similar probability: in fact the SN la 
rate in a galaxy ~ 13 Gyr after the beginning of star formation is 
given by: 

h la {\3)=k u xL c ^ AWdr + ^pjT^ JWr)drJ (4) 

where ki a is the number of SN la from a stellar population of unitary 
mass (or SN la productivity), and if/ c and if/ p are the (/i a weighted) 
average SF rates over, respectively, the last 0. 1 Gyr and from 0. 1 up 
to 13 Gyr ago. Equation l[4j explicitly shows the two contributions 
of prompt and tardy components of the DTD: if their ratio is of 0.3 
to 0.7, in order to have an equal probability to realise SD or DD 
explosions it suffice that the current SF rate exceeds the average SF 
rate in the past by a factor of ~ 2. 

To summarise, both the Solomonic and the Segregated mix- 
tures provide SD and DD explosions in virtually all galaxy types, 
but in the very old galaxies, where (almost) all events should come 
from the DD channel. It is worth noting that this conclusion only 
applies to Chandra explosions, since Sub-Chandra events could be 
common at very late delay times both in the Single and Double 
Degenerate scenarios. 

The scaling of the SN Ia rate in different stellar systems of- 
fers the opportunity of constraining the DTD, hence the SN Ia pro- 



genitor model (e.g. Paper i lMannucci et al .120061 ; IB! anc & G reggio 
2008). To this end I consider two sets of data: (i) the relation be- 
tween the SN Ia rate per unit luminosity and the mo rphological 
type of the parent galaxy (Greggio & Cappellaro 2009, hereinafter 
GC09) and (ii) the relation between the SN Ia rate per unit mass 
and th e specific SF rate (SSFR) of the parent galaxy dSullivan et al.l 
2006). Both relations (i) and (ii) basically stem from DTD being 
more populated at early delay times: as a consequence, the rate per 
unit mass is higher the younger the stellar system. A later galaxy 
type and/or a higher specific SF rate correspond to a larger fraction 
of young stars, hence a higher SN Ia rate per unit mass, or per unit 
luminosity, as far as luminosity traces the stellar mass. Nonetheless 
the two observational relations are largely independent (i) because 
of the different galaxy sample, and (ii) because of the different nor- 
malisation of the SN Ia rate, and different galaxy age tracer. There- 
fore, they yield complementary information on the DTD. In the 
next sections, the quantitative fit of the two relations is attempted 
separately. 



4.1 Rates as function of galaxy type 

The indication that late type gal axies are more effi cient in produc- 
ing SN Ia events goes back to iTammannl q 19781) . and was inter- 
preted a s due to a DTD more po pula ted at short delay times al - 
ready bv lOemler & Tinslevl i 19791) and lGreggio & Renzinil (l983). 
Since then, this observational feature has been put on firm grounds 
as a correlation between the SN Ia rate per unit B-band luminos- 
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Figure 10. Calibration of the parameter tsf of the SF rat e versus morpho- 
logica l galaxy type. The red dots show the galaxies of the Cappell aro et alj 
11999ft sample; the letters show the average location in color and type 
of the labelled galaxy class. The lines connect galaxy models obtained 
convolving the SF rate given by Eq. j5J with the tsf value labelled 
on the top axis (in Gyr) with simple stellar population (SSP) models. 
The latter were derived using the CMD tool (version 2.2) posted on 
the web (http://stev.oapd.inaf.it/cgi-bin/cmd i, with the followin g input: 
Marigo et al. 1 20()|) set of tracks; iMai'z Appellanizl feOOd) and iBessehl 
11990) photometric system; no circumstellar dust and no extinction; 
Salpeter IMF. The line type encodes the metallicity of the SSP models, as 
labelled. The blue crosses shows the colors of a model with an exponen- 
tially increasing SF rate in the last 13 Gyr, and an e-folding time of 1 Gyr, 
assumed to represent a bursting galaxy. 



ity (rate in SNu's Q) with the U - V colour of the parent galaxy 
( Cappellar o. Evans & Turattdl 19991) , and a correlation between the 
rate per unit K band lumino sity (rate in SNuK's) and the B - K 
colour of the parent galaxy (Mannucci et al. 2005). Both correla- 
tions are based on the (Cappellaro et al Jl 19991) sample of (nearby) 
galaxies. Ilvlannucci et aL 1 20051) also derived the rate per unit mass 
(rate in SNuM's) as a funct ion of the parent gal axy B - K colour, 
later used in Paper I and in Mann ucci et all |2006) to derive clues on 
the characteristics of the DTD. In particular, Paper I found that both 
CLOSE and WIDE DD models could explain the data, provided 
that the DTD slope was neither very flat nor very steep. However, 
the rate in SNuM's relies on theoretical models which are needed to 
convert the galaxies K-band luminosities into their mass; therefore 
the purely observed relations involve the rates per uni t lumi nosi- 
ties. The attempt to ac count for both ICappellaro et al. I Jl999l) and 
Mann ucci et all |2005) trends of the rates in SNu's and SNuK's 
with the parent galaxy U — V and B - K colors fails (GC09), with 
the rate per unit Lb requiring a natter DTD than the rate per unit 
L K . This is because the slope of the observed correlations depend 
on both the DTD and on the SF history assumed to model the galax- 
ies. To improve on the description of the latter, GC09 constrained 



Supernovae rate units are SNu's, SNuK's and SNuM's, corresponding to 
events per century per 10 10 Lb,o. Lk,o and M G respectively. 



the sample galax ies to be represented by the family of SF history 
laws adopted bv lGavazzi etail d2002l) : 



t / -f 

m = - r ex P \— r 

* cc \ * ox 



(5) 



which peaks at progressively older ages as the parameter r SF de- 
creases. A short tsf also implies a narrow age distribution, akin 
to early type galaxies, while for tsf > 13 Gyr the SF rate mono- 
tonically increases up to the present, having adopted an age of 13 
Gyr for all galaxies. A calibration of tsf vers us the morphological 
galaxy type T (as given in the RC3 catalogue. lde Vaucouleurs et all 
1991) was obtained by fitting the U — V and B - K colours of 
the sample galaxies. In this way theoretical SN la rates (in SNu's 
and in SNuK's) as a function of T could be computed and com- 
pared to the measured correlations. Such procedure has two advan- 
tages: the galaxies' SF history is constrained from their luminosity 
in four photometric ba nds, and the full sample of SN la events in 
Cappellaro et al. \ 19991) can be used to derive the rates in SNu's 
because all galaxies have a determined Lb and type T, while Lk is 
available for only a fraction of them. 

Fig. [10] shows the calibration of the tsf parameter versus T: 
the average colors of E, SO, Sa, Sb and Sc galaxies are reproduced 
with respectively r SF - 1,3. 4.5, 5.5 and 6 Gyr. It also appears that 
late type galaxies are better reproduced with a lower metallicity 
than early types. Concerning Irregular galaxies, the models account 
well for their average U—V color, but appear too red in B - K com- 
pared to the data. This may result from an inadequate description 
of their SF history, which could be dominated by a current burst. 
The cross in Fig.UOIshows the location of a bursting model, which 
well reproduces the B - K color of this galaxy type, but turns out 
too blue in U - V, hence not providing a satisfactory fit either. For 
all other galaxy types, however, the average colors are well repro- 
duced by the adopted SF history law, which is used in the following 
to fit the correlation of the SNIa rate with T. 

The black dots in Fig. II II shows the SNIa rates as function of 
morphological type determined in GC09, while the star shows the 
rate de termined for the merg ed class of early type galaxies (E and 
SO) in ICappellaro et al. ( 1999) (upper panel) and M annucci et all 
(lower panel). The latter value turns out higher than the rates 
determined by GC09 in the individual classes; this is mostly due to 
a (slightly) different binning of the parameter T which identifies the 
morphological types. Coloured lines in Fig.QT]show theoretical re- 
lations obtained with selected DTD models (see caption), meant 
to show the dependence on the main parameters. The models are 
computed adopting the correspondence between tsf and T just dis- 
cussed, plus a relation between t sf and metallicity Z, with Z de- 
creasing from 0.008 to 0.001 as t$f goes from 1 to 20. The crosses 
in Fig. [TJJ show the location of the bursting model (cross in Fig. 

rro>. 

It is worth emphasising that, rather then simulating a large 
population of individual objects, the approach here aims at describ- 
ing the photometric properties and the SNIa rate of the average 
galaxy within each bin of T; the morphological type is assumed 
to trace the age mixture in the galaxies, and describe its changes 
across the Hubble sequence. 

The theoretical rates in Fig. [TJJ are computed as the product 
of the rate per unit mass and the mass-to-light ratio of the galaxy 
models: 



»ia(0 X' 3 ^ (13 - T )/ia<T)dr 



J X13-T)£ Bm (T)dr 



(6) 
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Figure 11. SN la rate per unit B (up per panels) and per unit K (lower panels) luminosity, as function of the morphological type of the parent galaxy. Black dots 
show the observed rates from Greggio & Cappellaro 1 2009); the stars shows the published rates for the E+S0 class; coloured lines show models calculated with 
different DTDs and the SF rate given by Eq.{5). Left panels: theoretical relations obtained with single channel DTDs. The SD Chandra (SD-Ch,m green line) 
model assumes a Salpeter distribution for mi , a flat distribution of the mass ratios, and an efficiency e = 1 . Blue lines are obtained with DD-CLOSE models; 
DDC1 and DDC3 assume a distribution of separations skewed at small values, and differ for the minimum m%, respectively of 2.5 and 2 M Q ; DDC2 assumes a 
flat distribution of the separations and a minimum mi of 2.5 M Q . Red lines are obtained with DD-WIDE models, DDW1 and DDW2 being computed with the 
same choice of parameters as DDC1 and DDC2. Central and right panels: correlations obtained with mixed DTDs in the Solomonic and Segregated flavours; 
colours and line types reflect the specific DD model used to construct the mixture with the encoding labelled in the top left panel. The coloured crosses show 
the rates for the bursting model in Fig. |l()| put at an arbitrary value of T. 



where £b(k) are the light-to-mass ratios of simple stellar popula- 
tions, either in the B or in the K band, as functions of age, and M SF 
is the total mass transformed into stars (i.e. integral of tft(t)) up to 
the current galaxy age, here assumed of 13 Gyr. When going from 
early to late type galaxies, the increase of the first factor (rate per 
unit mass) competes with the decrease of the second factor (mass 
to light ratio), which is more pronounced in the B, rather than K, 
band. This explains the large variation of the model rate in SNuK's 
compared to that in SNu's; it also explains why the rate in SNu's 
can even be predicted higher in early than in late type galaxies (see 
e.g. model DDW2). 

The models are normalised to the observed value at T = 3, in- 
termediate between Sa and Sb galaxies; the normalisation yields a 
value for fc Ia , i.e. the productivity of SN la from one unit mass stel- 
lar population, which depends on the assumed IMF. The observed 
rates imply fci a 10~ 3 M G ' with relatively little dependence on the 
DTD or the photometric band, for the Salpeter IMF adopted in the 
CMD tool (i.e. <p(m) cc m' 23S down to 0.01 M Q ). A bottom light 



IMF, with <f>{m) cc m~ 13 between 0.1 and 0.5 M Q and a Salpeter 
slope above 0.5 M Q , would yield fci a values larger by a factor of 
~ 3. A detailed discussion of the uncertainty on k Ja requires a thor- 
ough evaluation of the quality of the fits, which goes beyond the 
scope of this paper. It is however interesting to notice that this value 
of the SNIa productivity, compar es well with th e total Fe content 
of galaxy clusters. According to iRenzinil |2004) the Fe Mass-to- 
Light ratio of clusters is M Fc /L B ~ 0.015 in solar units, assuming a 
stellar mass-to-light ratio of M,/L B ~ 3.5. This value is appropri- 
ate for a bottom light IMF, for which the current mass in stars (at 
old ages) is ~ 0.6 of the total mass transformed into stars. There- 
fore, the Fe content in galaxy clusters requires an Fe productivity of 
M Fc /M SF = 0.015/(3.5/0.6) = 0.0026. Adopting fc Ia = 0.003M G ' 
and an Fe production of 0.6 M G per SNIa event, the Fe productiv- 
ity from SNIa's is of 0.0018, which accounts for ~ 70 % of the Fe 
content in galaxy clusters; the remaining 30 % would be produced 
by CC SNe. 
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Fig. [TT] shows that the models with a flat distribution of the 
separations (dotted lines) tend to over predict the rate in Ellipticals, 
in both DD CLOSE and DD WIDE cases. This problem is less- 
ened for the Solomonic mixture, due to the dramatic drop of the 
SD channel at long delay times. However, when considering the 
merged E+SO class, the DD CLOSE model with a flat distribution 
of A becomes compatible with the data. It is then very important 
to derive a robust measurement of the SN la rate separately for the 
two morphological types: in early type galaxies, a difference in av- 
erage age of ~2 Gyr has a strong effect on the SN la rate per unit 
mass, if a sizable contribution from the SD channel exists. The fact 
that the rate in SNu's in SO galaxies appears so much higher than in 
Es favours this interpretation, although the statistics is still too poor 
to draw such conclusion. In this respect, the SD model appears to 
best represent the data. However, the green lines in Fig.QT]are ob- 
tained with a choice for the accretion efficiency and distribution of 
the mass ratios which maximize the fraction of systems with long 
delay times; any other choice would lead to an underproduction of 
SNIa events in E galaxies. 

Considering now the rates in late type galaxies, all models are 
consistent with the data within the very large error bar, although 
the predicted rates are lower than the average observed values. The 
parameter T n-X which has some impact on the fraction of events at 
short delay times (see Sect. 3.2) does not appear to play a role in the 
fit (e.g. DDC1 vs DDC3). In a star bursting regime the SN la rate 
gets substantially increased (crosses), especially when the DTD is 
particularly skewed towards the short delay times. As mentioned 
above, the galaxy modelling presented here does not reproduce sat- 
isfactorily the colors of Irregulars with either tested shapes for the 
SF rate. A more precise assessment of the recent SF in galaxies of 
the latest types is crucial to derive information on the DTD from 
their observed SN la rates. 

As a general conclusion, the trend of the SN la rate as function 
of the morphological type suggests that the DTD ought to be signif- 
icantly decreasing at the very long delay times. This disfavours flat 
distributions of the separations in the case of DD progenitors. Pro- 
vided that, the fit obtained with the mixed models is not particularly 
better than that with single channel models. 

4.2 Rates as function of specific star formation 

ISullivan et all (2006) measured the SN la rate per unit mass as a 
function of the current star formation rate per uni t mass for the 
Super nova Legacy Survey (SNLS) galaxy sample dSullivan et al.1 
2005). The correlation is based on 124 SN la at redshifts between 
0.2 and 0.75. It is worth recalling that in this work the galaxy 
masses were determined via s pectral energy distribution (SEP ) 
fitting to PEGASE 2.2 models dFioc & Rocca-Volmerangelll997h 
with IKroupal J200lh IMF, and that the SSFR is measured as the ra- 
tio between the mean SF rate over the last 0.5 Gyr and the current 
stellar mass of the galaxy, resulting from the SED fitting. 

Black dots in Fig. 1 121 show the obse rved relation for the star 
forming galaxies in lSullivan et alj ( 120061) : coloured dots and lines 
are models constructed for a wide variety of SF histories, so as to 
illustrate the impact on the theoretical relation of the adopted SF 
law, and of the galaxy age. Given their redshift range, the SN la 
host galaxies should have current ages between ~ 7 and 1 1 Gyr, or 
younger. Technically, the theoretical specific SN la and SF rates are 
computed as: 



M 



: 210~ 



(8) 



M SF 
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W ~ T)/laCT)dT 
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where t is the current age and Msf is the mass transformed into 
stars up to now. Because of the mass return from dying stars, M SF 
is always higher then the current stellar mass in a galaxy, by a fac- 
tor of ~ 1.6, 1.3 for old and young stellar populations, respectively. 
These values hold for a bottom light IMF, such as Kroupa's, and 
for galaxy ages older than ~ 1 Gyr. Therefore, accounting for the 
mass return corresponds to shifting the models up and right by 
~ 0.1 — 0.2 dex, the correction being however different for each 
model, since it depends on age and SF history. Given the differ- 
ences in the galaxy modelling between Sulliv an et al. I J2006I) and 
this paper, ra ther than correcting each model, th e data points (read 
off Fig. 3 in lPritchet, Howell & SullivarJl2008l) have been shifted 
by a variable amount (see caption). The applied correction is meant 
to represent the average ratio of the current stellar mass to the total 
mass gone into stars in the five bins. As usual, the productivity Aj a 
in Eq. Q is a free parameter; however the models in Figl 121 have 
been all computed with £ !a = 10 _2 8 Mg' irrespective of the DTD, a 
value which turns out to describe very well the SN la rate in the bin 
at SSFR=10~ 10 yrs~'. A refined fit of the SN la productivity seems 
unworthy, since the models and the data are not completely com- 
patible as just discussed, and since the SF history of the galaxies is 
not well constrained. 

The slope of the theoretical correlation between the rate per 
unit mass and the SSFR appears virtually insensitive to the assumed 
SF history, for star forming galaxies, as already pointed out by 
Pritchet et ai] j2008l) . The models overlap each other along a locus 
which is well described by a power law along which the younger 
the galaxy (either because of its current or its average age), the 
higher its SSFR, and its SN la rate per unit mass. The exponent 
of the power law is sensitive to the model DTD; for the models 
plotted in Fig(T2]it ranges from 0.6 (SD model) to 0.8 (Segregated 
mixture with DDC1 model). In addition, the theoretical relation is 
found systematically steeper for the Segregated mixtures compared 
to Solomonic mixtures or Single channel models. 

In general, all DTDs yield and acceptable representation of the 
data. Taken at face value, the observed correlation seems to favour 
models with a moderate power at the early delay times, as the SD 
model, or the Solomonic mixtures. The Segregated mixtures tend 
to predict a too high SN la rate in star bursting galaxies at given 
SSFR v alue, but no strong co nclusion can be drawn from this com- 
parison. lSullivan et alJd2006l) also determine the SN la rate per unit 
mass in passive galaxies (from SNLS data), which turns out very 
close to the rate measured for galaxies in the bin at SSFR=10" 
M /yr. This point is not plotted on Fig[l2] where, instead, I show 
the location of two models for old stellar populations (see caption), 
in which most of the stellar mass has been made at high redshift. 
These models have virtually null current specific SF rate and yet 
their current SN la rate is relatively high due to the late delay times 
compone nt of the DTD . This interpretation differs from what con- 
cluded in lPritchet et al.l d2008l) mainly because of the different ren- 
dition of the SF history in passive galaxies. Adopting Eq.|5]with a 
short t sf , the SSFR quickly drops as the time from the beginning 
of star formation increases, so that negligible SSFR are obtained 
when the SNIa rate is still relatively high because of the tardy com- 
ponent. For example, a model with r SF = 2 Gyr, at an age of 10 Gyr 
has a SSFR of 1.8 10~ 14 M Q /yr and an average age of 7.5 Gyr. With 
this star formation history, the SD DTD plotted in Fig.[l2]predicts a 
SNIa rate of 4.7 1 0~ 14 (Af a yr)~', match ing the rate in passive galax- 
ies determined bv lSullivanetai1d2006l) . However, as noticed in the 
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Figur e 12. C o rrelati on between the SN la rate and the SF rate both normalised to the total mass transformed into star. The data (black dots) are from 
I Sulliva n et al. 12006) shifted down and left to account for the mass return, by (0.12,0.12,0.15,0.17,0.19) dex for the 5 bins from right to left. The models 
(coloured lines and dots) assume different DTDs as labelled in each panel (see caption to Fig ll ll for the meaning of the acronyms). The colours encode the 
kind of SF history: (i) currently bursting, with a ratio of current to past SF rate of 500, 50 (light green), 10, 5 (orange); (ii) power laws with exponent -0.9 and 
-0.5 (blue), 0.5 and 0.9 (cyan); (iii) if/ given by Eq. J5J with t$f = 1 (magenta), 3, 5, 6 and 20 Gyr (red). The lines show the models as they age from 1 to 13 
Gyr; crosses show the location at 2 and 5 Gyr, full dots at 6, 7 and 1 1 Gyr, the stars mark the position at 13 Gyr. The two open dots show the rates of the model 



with ijj given by Eq. ^5} and tsf = 1 Gyr at an age of 7 (upper) and 1 1 (lower) Gyr. Notice that the SSFR of these models is much lower than 10 
virtually null. 



yr 



',and 



previous section, such SD model assumes a 100 % accretion of the 
envelope of the secondary onto the CO WD, and a flat distribution 
of q, so that the rate at late epochs is enhanced. For this reason, DD 
models, and/or Solomonic mixtures could be more plausible. 



5 SUMMARY AND DISCUSSION 

This paper presents the theoretical expectations for SN la progeni- 
tors focusing on two questions: the shape of the distribution of the 
delay times and the possible diversity of SN la progenitors of events 
at different delay times. These two issues impact on the cosmolog- 
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log DT (yr) 

Figure 13. Cumulative fraction of events as the delay time grows for the 
DTDs used in Fig.[TT] Dashed red, solid black and solid blue lines refer re- 
spectively to single channel, Solomonic mixtures and Segregated mixtures. 

ical use of SN la's as distance indicators because they may entail a 
systematic diversity between events at low and high redshift, where 
the progenitors are younger and more massive. As in Paper I, the 
approach refrains from describing the evolution of individual sys- 
tems, rather concentrating on the clock of the explosion and on a 
few parameters singled out as crucial for the two issues. In particu- 
lar, different from BPS renditions, the mass exchange phases during 
the binary evolution are not followed, and their result in promoting 
or not the final event is considered as an equally plausible option. 
The dependence of the likelihood of the SN la event on the initial 
binary separation is also neglected, having assumed that (i) there 
is a range of separations suitable for the realisation of the SN la, 
and that (ii) this range does not depend on the other binary param- 
eters (mi and m-x). While strictly speaking the latter assumption is 
untrue (see discussion at the beginning of Section 2.3), the compar- 
ison of the analytic DTDs with the results of BPS simulations seem 
to indicate that the crucial parameters are indeed the masses of the 
components. 

I have illustrated the main features of the analytic DTD for 
the single and the double degenerate evolutionary channels, con- 
centrating on Chandrassekhar exploders, since Sub-Chandra mod- 
els are currently disfavoured. Nevertheless, it is worth pointing out 
that Sub-Chandra events could be important especially to account 
for SN la's in old stellar populations. Actually these events are an 
embarrassment, since they are much easier realised than Chandra 
explosions, yet there are no established observational counterparts. 
Crucial parameters for shaping the SD DTD are the limits on the 
primary and on the secondary mass, and the efficiency of accretion 
onto the CO WD. The DTD is populated only within delay times 
comprised between the MS lifetime of the most and least massive 
secondary in the progenitor's system; a lower accretion efficiency 
results into a steepening of the DTD at late epochs. If only half of 
the envelope of the secondary is used for the mass growth of the 
CO WD, the DTD drops by three order of magnitudes between 1 
and 10 Gyr, making the rate in old stellar systems negligibly small. 



For both flavours of DD models (CLOSE and WIDE), the major 
parameters are the minimum mass of the secondary in successful 
systems and the slope of the distribution of the separations of the 
DD systems at birth. Additional parameters are the minimum val- 
ues for the delays t„ and tqw, the former being directly connected 
to the maximum value of m 2 in successful systems. The compari- 
son of the analytic DTDs with results in the literature shows a very 
good agreement with BPS renditions for the DD CLOSE channel 
(renditions for the DD WIDE channel are lacking in the literature). 
Less favourable is the comparison for the SD channel, for which 
the BPS results yield a large variety of DTD shapes. Clearly, in the 
simulations, the efficiency with which this channel produces SN la 
events is very sensitive to the criteria adopted to describe the mass 
exchange phases and the fate of the accreting WDs. 

5.1 Prompt and tardy components 

In principle, all channels accommodate prompt and tardy events, 
but the DTD is very continuous so that this distinction appears arti- 
ficial. As illustrated in Fig.[T3] for the single channel DTDs, the cu- 
mulative fraction of explosions smoothly increases with the delay 
time, so that the fraction of prompt events depends on the definition 
of the timescale over which the events are considered as such. As 
shown in Figsll llandll2lthese DTDs, with no strong discontinuity, 
do account well for the observed trends of the SN la rates with the 
parent stellar population. 

Motivated by recent observations which suggest that both SD 
and DD explosions may be realised in nature, I have constructed 
mixed models under two extreme assumptions: the Solomonic op- 
tion, in which the two channels contribute the same fraction of 
events within a Hubble time, and the Segregated mix, in which the 
SD and DD channels feed distinct ranges of delay times. Since the 
SD model has difficulties to account for late epoch explosions, I as- 
sumed that it provides a very prompt component of the DTD. This 
composition allows one to introduce a sharp discontinuity when the 
channel producing the events changes; the position in delay time of 
this discontinuity has been (arbitrarily) fixed a t 0.1 Gyr to mimic 
the DTD proposed by iMannucci et alj J20061) . In the attempt of 
constraining the DTDs, either single channel or mixed models, I 
have considered the observed correlations between (i) the SN la 
rate per unit luminosity and the parent galaxy type (GC09), and (ii) 
the SN la rate per unit mass and the specific SF rate of the parent 
galaxy dSullivan et alj 2006). Theoretical expectations have been 
computed with a population synthesis technique, adopting a variety 
of SF histories for the galaxies. To construct the first of the two cor- 
relations, the galaxies' (U - V) and (B - K) colours have been used 
to constrain the SF history; correlation (ii) turns out largely inde- 
pende nt of the adopted shape for the SF, as found by Pritc het et alj 
(2008). The observed correlation between the SN la rate per unit 
luminosity and the parent galaxy morphological type appears to re- 
quire a substantial depletion of the DTD at late delay times. This 
holds especially for the rate measured in SNu's and derives from the 
drop of the rate going from SO to E galaxies, which is not detectable 
when the two classes are merged together. In fact, the colours of SOs 
and Es indicate that the former are slightly younger, although both 
galaxy types are dominated by old stellar populations. SOs and Es, 
then, probe the DTD at very late delay times, where, e.g., the SD 
model predicts a considerable drop. The steepening of the rate go- 
ing from SO to Es measured in GC09 supports DTDs with a steep 
decline at late delays, as the SD channel, or Solomonic mixtures. 
The DD channel and the Segregated mix could also meet this con- 
straint, provided that the distribution of the separations is skewed at 
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the short end. This interpretation needs to be confirmed with larger 
galaxy samples and more SN la events. 

The correlation between the SN la rate and the galaxy type 
does not provide a robust constraint on the shape of the DTD at 
short delay times, because the SF history in the latest galaxies is 
not known with sufficient accuracy, and, admittedly, the colours of 
Irregulars are not well reproduced by the populations synthesis pre- 
sented here. The theoretical rates are found consistent with the data 
for Irregulars and Sc galaxies even with the relatively low fraction 
of prompt events of the SD model, but a better understanding of 
the recent and past SF history in these galaxies is needed to derive 
robust conclusions. 

The correlation between the SN la rate per unit mass and 
the SSFR is compatible with the models for most of the explored 
DTDs, although the rates in the star bursting galaxies indicate that 
the DTD should not be too skewed at the short delay times, thus 
disfavouring very steep DD CLOSE models and Segregated mix- 
tures. Models computed with the SD DTD yield a very nice de- 
scription of the correlation for the star forming galaxies; the rate 
in passive galaxies may also be well reproduced with this DTD if 
their inhabiting stellar populations have a (small) age spread: e.g., 
one model galaxy having formed ~ 80% of its current mass within 
the first 1.7 Gyr, at an age of 9 Gyr has a SN la rate|f| of ^ 4.2 1(T 14 
(M e yr) -1 very close to what determined for the passive galaxies in 
ISullivanet all {2006) (~ 5 10 14 (M Q yr) '). However, this estimate 
assumes that all of the envelope of the secondary is accreted by the 
CO WD companion, which looks contrived. 

5.2 SN la productivity 

The fit of the observed relations yields a value for the productivity 
fci a , i.e. the number of SN la's produced by a unit mass stellar pop- 
ulation over the whole range of delay times. In general ki a depends 
on the IMF, DTD and SF history assumed to model the galaxies, 
but the dependence on the DTD is small for systems in which the 
SF rate varies little over the whole range of ages, and actually van- 
ishes for a constant SF rate (see Paper I). Therefore, in this paper, 
the value of £ Ia is derived by matching the level of the SN la rate 
in galaxies of intermediate type (Sa-Sb). Assuming a bottom light 
IMF (Salpeter flattened below 0.5 M G ), from the rates in SNu's 
(SNuK's) in GC09 I find k la = 3.3(2.8) x IO^Mq 1 , with a 50 % 
total variation for all the DTDs plotted in Fig.[TT](single and mixed 
channels). This value requires that = 8 % of all stars with mass 
between 2.5 and 8 M Q produce a SN la event. For the Segregated 
mixture, 30 % of the progenitors are assumed to come from sys- 
tems with m 2 between 5 and 8 M Q : this corresponds to a likelihood 
of the SN la event of ~ 10% from the more massive progenitors, 
and of ~ 4% for stars with mass between 5 and 2.5 M Q . 

It is instructive to compare the SN la to the core collapse (CC) 
SNe productivity (k C c) for the same IMF. The number of CC events 
per unit mass from a stellar generation depends on mass limits of 
the progenitors; assuming an upper limit of 40 M Q , it amounts to 
6.5 (9) events every 1000 M Q for a lower limit of 10 (8) M Q , having 
adopted a Salpeter flattened IMF. Thus, each stellar generation is 
expected to produce ~ 2-3 CC SNe per SN la; the events, though, 
occur at different delay times: all the CC explosions take place soon 
after the star formation episode, while the SN la events are diluted 
over the whole Hubble time, with only the prompt component ex- 
ploding early. For the Segregated mixture considered here, a 1000 

5 This estimate assumes k\ d = 10~ 2 8 Mq' and 37 % a mass return. 



M Q stellar populations provides ^ 1 (prompt) SN la within 0. 1 Gyr 
from birth; the same stellar population provides ~ 7 CC SNe within 
30 Myr from birth. Therefore, a sample of elliptical galaxies sup- 
posed to have experienced a starburst within the last 0. 1 Gyr should 
yield 7 CC to 1 la prompt even ts; thus the excess of 7 SN la's in 
radio loud ellipticals noticed bv lDella Vallee t al. (2005) should be 
accompanied by ~ 50 CC events in the same galaxy sample, while 
none was observed. This observational result would require a de- 
tection efficiency of CC events ~ 50 times lower than that of la 
events, which is very hard to understand. The problem of the miss- 
ing CC SNe in radio loud Ellipticals becomes worse if the CC pro- 
ductivity is higher, as implied, e.g., by a lower limit of 8 M Q for 
their progenitors. It is worth noticing that although dust obscura- 
tion may preferentially hide CC SNe, since they explode soon after 
the onset of SF, in late type spirals 2.5 CC events a re detected ev- 
ery 1 SNIa ( Cappellaro, Botticella & Greagio 2007), which shows 
that dust does not have that dramatic effect at least in this kind of 
galaxies. 

The productivity ki- d ca n also be determined by fitting the rates 
from ISullivan et al. I fe006l) . Matching the SN la rate in galaxies 
with SSFR = 10- 10 (M yr)-' requires k la 1.6 x lO^M" 1 , which 
is ~ 2 times lower than what estimated from the rates per unit lu- 
minosity. Notice that with this lower productivity, the problem of 
the missing CC events in radio loud Ellipticals worsens. The dis- 
crepancy between the productivities estimated from the two kinds 
of correlations may be due to the different astrophysi cal ingredients 
for the stellar population models adopted here and in Sulliv anet al.l 
(2006): the IMF, the SF history (i.e. age and metallicity distribu- 
tion) of the galaxies, the stellar tracks to describe their light. In 
fact, in order to obtain the observational SN la rates and SSFR, as- 
sumptions are necessary for all of these ingredients. In addition, the 
estimates for k^ derived here represent zero point determinations, 
rather than accurate fits to the data, and the true distance between 
the best fitting k\ a values from the two correlations could be smaller 
than what quoted above. All in all, the results on the productivity of 
SN la's from matching the observed rates are rather encouraging, 
though there seems to be a systematic difference between the values 
derived from the two data sets, with the rate per unit luminosity re- 
quiring a factor of ~ 2 higher productivity. Whether this is a robust 
conclusion or not can be established only with complete consis- 
tency of assumptions in modelling and observational estimates. 

5.3 Diversity of SN la events 

The astrophysical parameter responsible for the diversity of 56 Ni 
production, and then of SN la light curves, has not been identi- 
fied yet. A number of possibilities have been proposed in the lit- 
erature, some of which related to the progenitor, like metallicity, 
age, ignition density, others more rel ated to the physical phenom- 
ena taking place during the expl osion ( Hoefl ich & Khokhlovll 19961: 
iKasen. Roepke & Wooslev8 2009). From an empirical point of view, 
the most conspicuous evidence concerns the correlation between 
the SN la light curve properties and the (average) age of the par - 
ent stellar population dGallagher et alj|2008l : iHowell et ail [2009). 
metallicity playing, in case, a secondary role. This suggests the pos- 
sibility that the diversity of SN la's is strongly related to their delay 
time. 

According to stellar evolution one expects a trend for the av- 
erage mass and temperature of the WD progenitor with delay time, 
but with a substantial spread. For the SD model, the average mass of 
the CO WD at the start of accretion first decreases, reaches a min- 
imum, but then increases again, since at the late delay times only 
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Figure 14. Differential (upper panels) and cumulative (lower panels) distributions of the delay times of events exploding in 13 Gyr old galaxies. Different 
panels show the results for different DTDs, as labelled in the upper figure; the mixed models assume a DD component of the CLOSE variety with r„ x =l and 
/? e =-0.975 (DDC3). The colour encodes the SF history, given by Eq|5]with different values of Tsf, labelled in the upper left panel. The values 20, 6, 5, 4, 3, 
and 1 Gyr produce model galaxies with colours akin to Irregulars, Sc, Sb, Sa, SO and E, respectively. For the Segregated mixture, the fraction of events with 
delay time shorter than 0. 1 Gyr obtained with the different SF histories is labelled in the upper right panel. The curves in the lower left panel are labelled with 
the (B luminosity weighted) average age of the galaxy model. Notice that the delay time scale in the upper panels has been expanded between and 1 Gyr for 
better readability of the figure. 



massive CO WDs are able to grow up to the Chandrassekhar limit. 
The range encompassed by the mass of the CO WD is maximum at 
T ~ 1 Gyr, but becomes narrower at late delays. For the DD WIDE 
channel, the total mass of the successful DD system is expected 
to vary between 1.4 and 2.4 M B (almost) irrespective of the delay 
time, since for each combination of mi and mi, small and wide 
separations should be realised, thereby providing events with any 
gravitational wave radiation delay. For the DD CLOSE channel, in- 
stead, we expect a trend of decreasing average Mod as t increases, 
because the more massive the progenitors, the smaller the sepa- 
ration of the DD system, and then the shorter the time needed to 
merge. The time between the formation of the CO WD and its final 
explosion is an increasing function of the delay time. Thus, for the 



DD CLOSE channel, the progenitors of events at short delay times 
are massive and hot (when accretion starts), while at long delay 
times they are of low mass and cool. According to Lesaff re et al.l 
(2006) models these two parameters could induce compensating 
trends on the ignition density, and the 56 Ni mass synthesised in the 
explosion may not show big variations over the total range of delay 
times. It is worth pointing out, however, that the models explore a 
range of relatively short cooling times (up to 1 Gyr), while the WD 
progenitors of SN la's in early type galaxies have been cooling for 
many Gys. The outcome of the explosion of CO WDs which start 
accreting when their structure is so highly degenerate needs to be 
investigated. 

Regardless the key astrophysical variable which causes the di- 
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versity of SN la light curves (whether ignition conditions, metallic- 
ity, evolutionary channel leading to the explosion), if this variable 
scales systematically with the delay time one expects a trend of the 
typical properties of SN la light curve with the age distribution in 
the parent stellar population. The observed correlation will depend 
on the SF history of the galaxies and on the DTD: the younger 
the stellar populations and/or the higher the fraction of short delay 
times, the higher the likelihood that a detected SN la event has a 
short t. 

As an illustrative example, Fig[l4]shows the distribution of the 
delay times of events occurring in galaxies with 13 Gyr of age, hav- 
ing experienced different SF histories. The differential distributions 
are computed as: 

X r+A >(13-r)/ la (r)dr 
An = — (9) 

/ >(13-T)/ hl (T)dT 

and represent the probability that a SN la detected in a galaxy of a 
certain type has a given delay time. It appears that in the late galaxy 
types most of the events cover a range of delay times up to about 
1-2 Gyr: as shown in Sect. 3 this is wide enough to encompass 
the whole range of masses of the CO WDs. Only in the case of 
strongly star bursting galaxies the majority of events may come 
from a truly prompt component, i.e. the most massive and hottest 
CO WDs. Even for the Segregated mixture in combination with 
the model for Irregular galaxies, more than 40 % of the SN la's 
should have a delay time in excess of 0.5 Gyr. In Sa galaxies all the 
delay times have about the same probability of being sampled, and 
the distribution remains wide also in SO galaxies, though mostly 
confined at t > 5 Gyr; in Ellipticals, instead, we expect to sample 
only the latest delay times. Therefore, if diversity were related to 
the delay time, SN la's in Es should be rather uniform, while in 
later galaxy types we should detect a large variety of light curve 
properties, as suggested by the data. Notice that in the DD WIDE 
channel the mass range of the successful DD systems is not likely 
to become narrower at late delays: for this model, if diversity was 
connected to Mdd, a wide range of light curve properties should be 
observed also in the earliest galaxy types. 

Uniformity of events in old galaxies is supported in 
iHowell et alj {2009), but it is difficult to understand quantitatively 
the result that galaxies with a light averaged age in excess of 3 
Gyr host only faint supernovae (technically with a low ejected 56 Ni 
mass). The curves plotted in the bottom left panel of Fig[l4]are la- 
belled with the B-light averaged age of each galaxy model. Only 
Irregulars and Sc galaxies appear younger than 3 Gyr, while Sb's 
and Sa's have an average age of several Gyrs. In these galaxies 
the distribution of the delay times of their hosted SN la's is pre- 
dicted rather broad, which would imply a broad range of photomet- 
ric properties, irrespective of the specific DTD model. Therefore, 
it is difficult to reconcile the notion that delay time is primarily re - 
sponsible for diversity of SN la events with lHowell et alj (2009) 
result. 



6 CONCLUSIONS 

According to expectations of stellar evolution in close binaries the 
DTD is a continuous distribution, skewed at the short delay times, 
but with no real justification for a separate prompt component. 
Mixed models may implement this concept, e.g. considering that 
the two channels (SD and DD) feed distinct ranges of delay times. 
However, there's no clear theoretical reason why it should be so, 



and the observed correlations between the specific rates and the 
parent galaxy properties do not appear to require it. 

Indeed, the observed correlations do not favour mixed over 
single channel models. The indication of the two channels being 
active comes from the data of individual explosions. More data are 
needed to establish whether the SD explosions are favoured at short 
delay times, i.e. in young stellar populations, which could support 
the Segregated mixture models presented here. 

The correlation between the SN la rate per unit B-luminosity 
and the parent galaxy type suggests a steep drop of the DTD at late 
epochs, consistent with the SD or DD models with a distribution of 
separations depleted of wide systems. The correlation between the 
SN la rate per unit mass and the specific SF rate is (broadly) con- 
sistent with this kind of models, and seems to require a relatively 
small fraction of prompt events. However, in order to constrain the 
fraction of short delay times in the DTD we need a more accurate 
description of the SF history in galaxies with high SSFR and of the 
latest morphological types. 

The observed rates allow us to evaluate the SN la productivity 
£i a , which appears to be in the range (1.6-3.3) 10~ 3 Mq'. Although 
it is reassuring that the different data sets yield values which agree 
within a factor of ~ 2, there may be a systematic difference between 
the productivity evaluated from the rates per unit luminosity and 
per unit SSFR. On the other hand, the ingredients adopted for the 
galaxies' modelling to fit the two correlations are not exactly the 
same. Therefore the discrepancy on the SN la productivity may 
just reflect a theoretical uncertainty. 

Concerning the diversity of SN la light curves, we expect a 
wide range for the mass of the CO WDs at the start of the accretion 
phase in most galaxy types; only in the older stellar populations the 
properties of the progenitors should be more uniform. For mixed 
models, in both the Solomonic and Segregated flavours, late type 
galaxies should host SD and DD events with comparable probabil- 
ity, while in early type galaxies most, if not all, SN la should come 
from DD progenitors. Only in very young stellar populations, with 
very high specific SF, the SD events are likely to prevail. Current 
data on the correlation between the light curve properties and the 
host galaxy seem to yield contradictory indications, possibly be- 
cause of the use of different parameters to characterise the SN la 
light curve and the properties of the host galaxy. Progress requires 
big samples of events, as well as more accurate knowledge of the 
age distribution in SN la hosts; on the theoretical side more models 
are also needed, in particular exploring the fate of an accreting CO 
WD which has been cooling for several Gyrs, hence starting from 
highly degenerate conditions. 

A broader diversity of progenitor's properties is expected in 
local SN la samples compared to high redshift ones, since the for- 
mer encompass all possible delay times, while the latter necessarily 
lack events with delay times longer than the corresponding age of 
the universe. Nevertheless, a fair fraction of explosions in nearby 
star forming galaxies (spirals and irregulars) should have a short de- 
lay time, e.g. shorter than 1 Gyr. Therefore the empirical standard- 
isation of the light curve obtained from the SN la's in nearby star 
forming galaxies should be applicable also to high redshift events, 
thus supporting SN la's as viable standard candles all the way to 
cosmological distances. 



ACKNOWLEDGMENTS 

I am indebted to Alvio Renzini for helpful discussions and com- 
ments, and for careful reading of the manuscript. I also thank En- 



The rates of Type la Supernovae. II. 



21 



rico Cappellaro for discussions and for checking on the discrep- 
ancy between the SN la rate in early type galaxies determined by 
iMannucci et al ] J2005h and by GC09. Finally, I thank the anony- 
mous referee for productive comments. 



REFERENCES 

Abt H.A., 1983, ARA&A, 21, 343 
Altavilla G. et al., 2004, MNRAS, 349, 1344 
Arnett D., 1982, ApJ, 253, 785 

Belczynski K., Bulik T., Ruiter A.J., 2005, ApJ, 629, 915 
Bessell M.S., 1990, Publ. Astron. Soc. Pac, 102, 1181 
Blanc G., Greggio L., 2008, New Astron., 13, 606 
Cappellaro E., Evans R., Turatto M., 1999, A&A, 351, 459 
Cappellaro E., Botticella M.T., Greggio L., 2007, in Immler S., 

Weiler K., McCray R., eds, AIP Conf. Proc. 937, Supernova 

1987A:20 Years After: Supernovae and Gamma-Ray Bursters. 

Am. Inst. Phys., New York, p. 198 
Ciotti L., D'Ercole A., Renzini A., Pellegrini S., 1991, ApJ, 376, 

380 

Delia Valle M., Panagia N., Padovani P., Cappellaro E., Mannucci 

E, Turatto M., 2005, ApJ, 629,750 
de Vaucouleurs G, de Vaucouleurs A., Corwin H.G, Buta R.J., 

Paturel G, Fouque P. , 1991. Third Reference Catalogue of 

Bright Galaxies. Springer- Verlag, New York 
Di Stefano R, Primini FA., Liu J., Kong A., Patel B., 2009, 

preprint flarXiv:0909.2046t 
Ellis, R.S. et al., 2008, ApJ, 674, 51 
Fioc M., Rocca-Volmerange B., 1997, A&A, 326, 950 
Gallagher J.S., Garnavich P.M., Caldwell N., Kirshner R.P, Jha 

S.W., Li W., Ganeshalingam M., Filippenko A,.V., 2008, ApJ, 

685, 752 

Gavazzi G., Bonfanti C, Sanvito G, Boselli A., Scodeggio M., 
2002, ApJ, 576, 135 

Greggio L., 1996, in Kunth D., Guiderdoni B., Heydari-Malayeri 
M., Thuan T.X., eds. The Interplay Between Massive Star For- 
mation, the ISM, and Galaxy Evolution. Edition Frontieres, Gif- 
sur-Yvette, p. 98 

Greggio L„ 2005, A&A, 441, 1055 (Paper I) 

Greggio L., Cappellaro E., 2009, in Antonelli L.A., Brocato E., 
Limongi M., Menci N., Raimondo G., Tornambe A., eds, AIP 
Conf. Proc. 1111. Probing Stellar Populations Out to the Distant 
Universe. Am. Inst. Phys., New York, p.477 (GC09) 

Greggio L., Renzini A., 1983, A&A, 118, 217 

Greggio L., Renzini A., Daddi E„ 2008, MNRAS, 388, 829 

Hachisu I., Kato M., Nomoto K., 1996, ApJ, 470, LI 00 

Hachisu I., Kato M., Luna G.J.M., 2007, ApJ, 659, L153 

Hachisu I., Kato M., Nomoto K, 2008, ApJ, 683, L127 

Han Z., Podsiadlowski Ph., 2004, MNRAS, 350, 1301 

Hicken M. et al. , 2009, ApJ, 700, 331 

Hillebrandt W„ Niemeyer J.C., 2000, ARA&A, 38, 191 

Hoeflich P., Khokhlov A., 1996, ApJ, 457, 500 

Hoeflich P., Khokhlov A., Wheeler J.C., Phillips M.M., Suntzeff 
N.B., Hamuy M., 1996, ApJ, 472, L81 

Howell D.A. et al., 2006, Nature, 443, 308 

Howell D.A. et al., 2009, ApJ, 691, 661 

Ibenl.Jr, 1991, ApJS, 76, 555 

Iben I.Jr, Tutukov A.V., 1987, ApJ, 313, 727 

Iben IJr, Tutukov A.V., 1994, ApJ, 431, 264 

Lesaffre P., Han Z., Tout C.A., Podsiadlowski Ph., Martin R.G., 
2006, MNRAS, 368, 187 



Kalirai J.S., Hansen B.M.S., Kelson D.D., Reitzel D.B., Rich 

R.M., Richer H.B., 2008, ApJ, 676, 594 
Kasen D., 2010, ApJ, 708, 1025 

Kasen D., Roepke F.K., Woosley S.E., 2009, Nature, 460, 869 
Kenyon S.J., Livio M., Mikolajewska J., Tout C.A., 1993, ApJ, 
407, L81 

Kouwenhoven M.B.N., Brown A.G.A., Portegies Zwart S.F., 

Kaper L., 2007, A&A, 474, 77 
Kroupa P., 2001, MNRAS, 322, 231 
Mai'z Appellaniz J., 2006, AJ, 131, 1184 

Mannucci E, Delia Valle M., Panagia N., Cappellaro E., Cresci 

G, Maiolino R., Petrosian A., Turatto M., 2005, A&A, 433, 807 
Mannucci E, Delia Valle M., Panagia N., 2006, MNRAS, 370, 

773 

Matteucci M.F., Greggio L., 1986, A&A, 154, 279 

Matteucci M.F., Panagia N., Pipino A., Mannucci E, Recchi S., 

Delia Valle M., 2006, MNRAS, 372, 265 
Marigo P., Girardi L., Bressan A., Groenewegen M.A.T., Silva L., 

Granato G.L., 2008, A&A, 482, 883 
Mazzali PA., Nomoto K, Cappellaro E., Nakamura T, Umeda 

H, Iwamoto K, 2001, ApJ, 547, 988 

Meng X., Chen X., Han Z., 2009, MNRAS, 395, 2103 
Mereghetti S., Tiengo A., Esposito P., La Palombara N., Israel 

G.L., Stella L., 2009, Sci, 325, p. 1222 
Munari U., Renzini A., 1992, ApJ, 397, L87 
Napiwotzki R., et al., 2005,in Koester D., Moehler S, eds, ASP 

Conf. Series 334. 14th European Workshop on White Dwarfs. 

Astron. Soc. Pac, San Francisco, p.375 
Neill J.D. et al., 2009, ApJ, 707, 1449 

Nelemans G., Verbunt E, Yungelson L.R., Portegies-Zwart S.F., 

2000, A&A, 360, 1011 
Nobili S., Goobar A., 2008, A&A, 487, 19 
Nomoto K, Iben I.Jr, 1985, ApJ, 297, 531 

Nomoto K, Kamiya Y, Nakasato N., Hachisu I., Kato M., 2009, 
in Antonelli LA., Brocato E., Limongi M., Menci N., Raimondo 
G., Tornambe A., eds, AIP Conf. Proc. 1111. Probing Stellar 
Populations Out to the Distant Universe. Am. Inst. Phys., New 
York, p.267 

Nugent P., Baron E., Branch D., Fisher A., Hauschildt PH., 1997, 

ApJ, 485, 812 
Oemler A., Tinsley B.M. , 1979, AJ, 84, 985 
Parthasarathy M., Branch D., Jeffery D.J., Baron E., 2007, New 

Astron. Rev., 51, 524 
Patat F. et al. , 1997, Sci., 317, 924 
Perlmutter S. et al. , 1997, ApJ, 483, 565 
Perlmutter S. et al. , 1999, ApJ, 517, 565 
Phillips M.M., 1993, ApJ, 413, L105 

Piersanti L., Tornambe A., Straniero O., Dommguez I., 2009, in 
Antonelli L.A., Brocato E., Limongi M., Menci N., Raimondo 
G., Tornambe A., eds, AIP Conf. Proc. 1111. Probing Stellar 
Populations Out to the Distant Universe. Am. Inst. Phys., New 
York, p.259 

Pinsonneault M.H., Stanek K.Z., 2006, ApJ, 639, L67 
Pritchet C.J., Howell D.A., Sullivan M., 2006, ApJ, 683, L25 
Riess A.G. et al. , 1998, AJ, 116, 1009 

Renzini A., 2004, in Mulchaey J.S., Dressier A., Oemler A., 
eds. Clusters of Galaxies: Probes of Cosmological Structure and 
Galaxy Evolution. Cambridge Univ. Press, Cambridge, p. 260 

Renzini A., 2006, ARA&A, 44, 141 

Ruiter A.J., Belczynski K, Fryer C, 2009, ApJ, 699, 2026 

Salaris M., Garcfa-Berro E., Hernandez M., Isern J., Saumon D., 
2000, ApJ, 544, 1036 



22 Laura Greggio 



Simon J.D. et al., 2009, ApJ, 702, 1157 

Sullivan M. et al. , 2005, in Turatto M., Benetti S., Zampieri L., 
Shea W., eds, ASP Conf. Series 342. 1604-2004: Supernovae 
as Cosmological Lighthouses. Astron. Soc. Pac, San Francisco, 
p.466. 

Sullivan M. et al. , 2006, ApJ, 648, 868 

Tammann G.A., 1977, Mem. Soc. Astron. Ital., 49, 315 

Timmes F.X., Brown E.F., Truran J.W., 2003, ApJ, 590, L83 

Trimble V., 2008, The Observatory, 128, 286 

van den Heuvel E.P.J. , Bhattacharya D., Nomoto K., Rappaport 

S., 1992, A&A, 262, 97 
Yamanaka M., et al., 2009, ApJ, 707, LI 18 
Yoon S.-C, Podsiadlowski Ph., Rosswog S., 2007, MNRAS, 380, 

933 

Yungelson L.R., 2005,in Burderi L., Antonelli L.A., D' Antona F, 
Di Salvo T, Israel G.L., Piersanti L., Tornambe A., Straniero 
O., eds, AIP Conf. Proc. 797. Interacting Binaries: Accretion, 
Evolution and Outcomes. Am. Inst. Phys., New York, p.l 

Yungelson L.R., Livio M., 2000, ApJ, 528, 108 

Wang B., Meng X., Chen X., Han Z., 2009, MNRAS, 395, 847 

Webbink R.F, 1984, ApJ, 227, 355 

Williams K.A., Bolte M., Koester D., 2004, ApJ, 615, L49 
Woosley S.E., Weaver T.A., 1994, ApJ, 423, 371 



